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Abstract. A general class of (multiple) hypergeometric type integrals asso- 
ciated with the Jacobi theta functions is defined. These integrals are related 
to theta hypergeometric series via the residue calculus. In the one variable 
case, theta function extensions of the Meijer function are obtained. A number 
of multiple generalizations of the elliptic beta integral IS2I associated with the 
root systems A n and C n is described. Some of the C n -examples were proposed 
earlier by van Diejen and the author, but other integrals are new. An example 
of the biorthogonality relations associated with the elliptic beta integrals is 
considered in detail. 
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1. Introduction 

Exact integration formulas and integral representations of functions are impor- 
tant from various points of view. Such representations serve sometimes as defini- 
tions of functions, but more often they are needed for the better understanding 
of properties of functions defined beforehand. Due to numerous applications (see 
AAR ), the Euler beta integral 



(l-xf^dx^ 



Re a, Re (3 > 0, 



(1.1) 



where T(z) is the standard gamma function, plays a fundamental role in classical 
analysis. Various ^-generalizations of (|1.1|) involving g-gamma functions have been 
proposed within the theory of basic hypergeometric series [GTT] . Recently, a "third 
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floor" of the hierarchy of beta-type integrals, which is related to the elliptic gamma 
function, has been discovered in |S2| (in the one variable case) and in |DS1I IDS2] 
(multiple extensions). For a brief review of results in this direction, see DS4 . 
In this paper, we discuss a general class of hypergeometric type integrals associ- 
ated with the Jacobi theta functions and propose several new multiple elliptic beta 
integrals admitting exact evaluations. 

An infinite hierarchy of multiple gamma functions was proposed by Barnes long 
time ago |Bal| : 



where 7^ are some constants analogous to the Euler constant and £1 — n\U}\ + 
. . . + n r uj r (if some of the ratios LUi/u>k are real, then they must be positive). The 
prime in the product sign means that the point rii = . . . = n r = is skipped. The 
function i|1.2fl satisfies a collection of r first order difference equations 

T r (lt + U>j] u>) 1 , 

T r (u;u) r r _i(u;o;(j))' 

where tu(j) = (u>i, . . . , Uj—i,u)j+i, . . . , u> r ) and u>) — p{uj)uj u / u) T(u/lo) for some 
constant p(ui) (for a brief account of this function, see also appendix A in [JM ) . 

Following Barnes' analysis, in jj] Jackson has considered the generalized gamma 
functions in a slightly different way and proposed the g-gamma function and the 
elliptic gamma function. We recall the definition of the latter. Taking two complex 
variables q and p such that \q\, \p\ < 1, we compose the following (convergent) 
Jackson double infinite product: 

oo 

(z\q,p)oo= J] (1-3) 

j,k=0 

Two first order q- and p-difference equations for this product 

- {z;p)oo, t ^ = (2;g)oo, (1.4) 



{qz'iq,p)oo ' (pz;q,p) 

[oo 
fe=0 



where (z; p)^ — Il^Lo(^ — z P h )i are °f major importance. Replacing z by pz~ 



the first equation and by qz in the second, we get 

— — = {pz ;p)oo, — = {qz ;q) 00 . (1.5) 

[qpz l ;q,p)oo (qpz l ;q,p)oo 

We define a theta function as follows: 

0(z;p) = (z;p) OQ (pz~ 1 ;p) oc . (1.6) 
It is related to the standard Jacobi (9i-function |WW| in a simple way 



oo 

0i(u;<r,T)=-i (-1)'V 2 " +1)V V" +1/2)U 

71 — — OO 

= p 1 / 8 iq- u / 2 (p;p) OD e(q u ;p), u G C, 

where we assume that g = e 27rlCT ,p = e 27r4r . Sometimes, for brevity, it is convenient 
to drop q and p or the modular parameters a and r in the notations for theta 
functions and elliptic gamma functions, as well as for the elliptic analogs of shifted 
factorials to be defined below. 
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The function 9i(u) is entire, odd, 6\(— u) = —6\{u), and doubly quasiperiodic 

6i(u + a' 1 ) = -0i (it), 0i (u + rcr- 1 ) = -e^- 2 **™^). (1.7) 

These transformation properties of the 0i -function are extensively used in our for- 
malism. For the 6(z;p) function, they take the form 

9(pz;p) = eiz-'-p) = -z-^iz-p). (1.8) 

Now we multiply the left-hand sides and right-hand sides of the first identities 
in (|1.4fl and (|1.5fl . respectively, and do the same with the second identities. This 
yields the difference equations 

F(qz; q,p) = 6(z;p)T{z\ q,p), T(pz; q,p) = 9(z; q)T(z; q,p), (1.9) 

for the elliptic gamma function T(z; q-,p); in the explicit form, we have 

T(z;q,p)= J] , \l ■ (1-10) 

j,fc=0 

Despite of the fact that the general idea of associating a generalized gamma function 
with the elliptic theta function was formulated in the well-known paper p] , it did not 
get much attention. However, Jackson's double infinite product was used explicitly 
in the mathematical physics literature on integrable models of statistical mechanics 
starting with the Baxter's work on the eight vertex model; see |Bax| . The name 
"elliptic gamma function" for the product (|1.10|) was proposed by Ruijsenaars in 
the recent paper |Rul| . where he reintroduced the function T(z;q,p) anew and 
started a systematic investigation of its properties. A further detailed analysis of 
this function was performed by Felder and Varchenko in |FV| . 

In order to compare the elliptic gamma function with the Barnes multiple gamma 
function, in l|1.10|l we put 

27ri— 27ri— 2iri^2- 

z = e q — e p — e ""2 , 

where u>i are some constants satisfying the same constraints for the quasiperiods as 
in ljl.2|l . Then it is not difficult to see that the set of zeros and poles of T(z; q,p), 
viewed as a meromorphic function of the variable u, coincides with the set of zeros 
and poles of the following combination of Barnes Ta-functions: 

r 3 (wi +W3 - u;cji, u>2, u>3)r 3 (o;i + lo 3 - lo 2 - u;lji, -lo 2 ,lu 3 ) 

This means that the ratio of q,p) and (|1.11|) is an entire function of u, and this 
function is seen to be given by an exponential of some polynomial of u of the third 
degree. 

For arbitrary complex s, the elliptic shifted factorials are defined as ratios of 
elliptic gamma functions 

Q , \ ^{zq s ;q,p) 

0{z;p;q) s = — r> 
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Ms(t)= - r 1 - <s - 4 4 rr. " ^ 



We use also the following shorthand notation: 

k 

T(t 1 ,...,t k ;q,p) = Y[T(tj;q,p), 

i=i 

n—1 

0(t 1 ,...,t k ;nq) n = Y[Y[0(t ]q e ;p) 7 n e N. 

0=1 e=o 

The elliptic beta integral |S2| is the first exact integration formula involving the 
elliptic gamma function. We conclude this section by an explicit description of it. 

Theorem 1. Let five complex parameters t m , m = 0, . . . ,4, satisfy the inequalities 
\t m \ < 1, \pq\ < \A\, where A = Ilr=o^ r - Define the elliptic beta integral as the 
following contour integral: 

f dz 

M E (t)= / A E (z,t) — , (1.12) 
Jt z 

where T is the positively oriented unit circle, and 

A£(Z,t) - 2^T{z\z-\zA lZ -^A-q,p)- (L13) 

Then 

2 llo< m < s <4 r (^^;g,rt 

<7)oo(p;p)oo rim=o r^™ 1 ; ?>p) 

Relation (|1.12|) determines a new Askey- Wilson type integral representing an 
elliptic extension of the Nassrallah-Rahman integral. Indeed, if we set p = 0, then 
the integral p. 12(1 is reduced to the Nassrallah-Rahman g-beta integral |NRI IR1| 
which, in turn, is a one parameter extension of the celebrated Askey- Wilson g-beta 
integral |AW| . Theorem 1 was proved by the author with the help of an elliptic 
generalization of the method used by Askey in |As| for proving the Nassrallah- 
Rahman integral. A large list of known plain and g-hypergeometric beta integrals 
was given in |RS2| . 

As shown in |DS1| . a special finite-dimensional reduction of (|1.12|l associated 
with the residue calculus results in the elliptic generalization of the Jackson sum 
for a terminating g&j basic hypergeometric series, which was discovered by Frenkel 
and Turaev in JFTJ . In |S4j , the integral H1.12|l was applied to the construction of a 
large family of continuous biorthogonal functions generalizing the Rahman's 10^9 
biorthogonal rational functions |R1I IR2| . These functions are expressed through 
products of two 12-E11 elliptic hypergeometric series with different modular param- 
eters (for the definition of an appropriate system of notations for such series, see 
|S5p. It is believed that in the theory of biorthogonal functions they play a role 
similar to that played by the Askey- Wilson polynomials |AW| in the theory of or- 
thogonal polynomials. We describe these biorthogonal functions in the last section 
and in Appendix A of the present paper and give complete proofs of some results 
announced in S 1 . An elliptic extension of Wilson's discrete (finite-dimensional) 
set of biorthogonal rational functions |Wij was constructed earlier by Zhedanov and 
the author in jSZll ISZ2) . The corresponding three terms recurrence relation gener- 
ates the most general example in the pool of known terminating continued fractions 
expressed in terms of the series of hypergeometric type, namely, it is expressed via 
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the 12-E11 series as well [SZ3I ISZ4] . The integral (|1.12|) leads to integral representa- 
tions for a terminating 12-E11 series and its particular bilinear form |S4| (the proofs 
are given in Appendix B). All these results open new ways of exploration of the 
world of elliptic functions and modular forms, which complement recent progress 
reached in the classical setting by Milne |Mi2j . 



Acknowledgments 

This paper is dedicated to Mizan Rahman, a master of g-special functions from 
whom the author has learned much over ten years of regular contacts. The author 
is also grateful to J.F. van Diejen for many useful discussions during the work on 
|DSl| - |Dlv4] . and for remarks to this paper. Permanent encouragement from A.S. 
Zhedanov is highly appreciated as well. Special thanks go to C. Krattenthaler for 
explaining some of the Warnaar theorems jWaj , to V. Tarasov for emphasizing the 
importance of the Barnes work |Bal[ and for drawing attention to the paper |TVj . 
and to S. Kharchev for a discussion of the properties of the generalized gamma 
functions. Some of the main results were obtained during the author's stay at the 
Max-Planck-Institut fiir Mathcmatik (Bonn) in the summer of 2002. The author is 
indebted to this institute for warm hospitality and to Yu.I. Manin, A. Okounkov, 
and D. Zagier for several stimulating discussions at MPI of the subject of the 
present paper. The organizers of the workshops "Special Functions in the Digital 
Age" (Minneapolis, USA, July 22-August 2, 2002) and "Classical and Quantum 
Integrable Systems" (Protvino, Russia, January 9-11, 2003) are thanked for giving 
an opportunity to present there some results of this work. This research is supported 
in part by the Russian Foundation for Basic Research (RFBR) Grant No. 03-01- 
00781. 

2. A GENERAL DEFINITION OF THETA HYPERGEOMETRIC INTEGRALS 

The right-hand side of (|1.12|) belongs to a general class of integrals related to 
the series of hypergeometric type built from Jacobi theta functions. In accordance 
with the theory of general theta hypergeometric series developed in |S5j . we give 
the following definition. 

Definition. Let C denote a smooth Jordan curve on the complex plane. Let 
A (3/1, . . . , y n ) be a meromorphic function of its arguments 3/1, . . . , y n . Consider the 
(multiple) integrals 



In = dyi... dy n A(y x , . . . , y n ) (2.1) 
Jc Jc 

and the ratios 

My) = A( A y ;'--" + 1 '--^ . (2.2) 
A (3/1, ...,ye,...,y n ) 

Then the integrals I n are called: 

1) the plain hypergeometric integrals if 

hi(y) = Ri(y), (2.3) 

are rational functions of 3/1, . . . , y n for all I = 1, . . . , n; 

2) the q-hypergeometric integrals if 

h t (y) = R e (g y ), (2.4) 
are rational functions of q yi , . . . , q Vn , q G. C, for all I = 1, . . . , n; 
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3) the elliptic hypergeometric integrals if for all I = 1, . . . , n the ratios he(y) are 
elliptic functions of the variables y%, . . . , y n with periods cr _1 and ra _1 , Im(r) > 0; 

4) the general theta hypergeometric integrals if he(y) and l/he(y) are meromor- 
phic functions obeying the double quasiperiodicity conditions 

ht(y u ...,y k + o-\ ...,y n ) = e^U ^Ww+'afc^), 

fcfc/!, . . . , y k + To-~\ ...,y n ) = e^U ^)vi+^ hi (y), (2.5) 

with the quasiperiodicity factors similar to those for the Weierstrass sigma function 
(which is related to 6\(u) in a simple way, see |\\ \\ j i. 

If we assume that the variables y%, . . . , y n are discrete, y S N", and replace inte- 
grals by sums X) y eN™ > then we get the definitions of the plain and g-hypergeometric 
series, which go back to Horn |OGR| . and the definition of the elliptic hypergeo- 
metric series suggested in |S5j . respectively. The theta hypergeometric series were 
defined in |S5j in a less general form because of the less general choice of quasiperi- 
odicity factors. Evidently, if aik(j) = bo- = cik{j) = dik = 0, then the theta 
hypergeometric functions are reduced to the elliptic ones. The integrals (or series) 
defined in this way do not form an algebra because, in general, sums of hypergeo- 
metric integrals do not fit the taken definition. 

The shifts yi — ► ye + 1 in (|2.2[l may be replaced by translations by an arbitrary 
constant yg — > yi + u)\, uj\ — const. However, we can replace oj\ by 1 after an 
appropriate rescaling of yg, which results in a simple deformation of the contour C 
in (|2~T)l . 

Consider the case of n = 1 in detail. The general rational function of y can be 
represented in the form 

( , _ IXU (1 - Qj + v) n^n+i ( a j - i - y) 

where n, r, m, s are arbitrary integers, x is an arbitrary complex constant, and aj, bj 
describe the positions of the zeros and poles of R(y). The equation A(y + 1) = 
R(y)A(y) has the following general solution: 

r(6,- - y) IT?-i T(l - aj + y) 
Ilj= m+ i r(l - bj + y) n j= „+i r(aj - y) 

where T(y) is the standard gamma function and f(y) is an arbitrary periodic func- 
tion, (p(y + 1) = ip(y). If we set f(y) — 1, then for an appropriate choice of the 
contour C the integral I\ (see H2.ll ) is none other than the Meijer function [EMOT . 
In this case we have no natural additional tools for fixing an infinite dimensional 
(functional) freedom contained in the solution A(y). 
In the g-case, in a similar way we can write 

[Q ) lir : : ' "'A ,,,.-1 -,'/")•''• 

For < \q\ < 1, the general meromorphic solution of the equation A(y + 1) = 
R(qy)A(y) is 
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where, again, tp(y) is an arbitrary periodic function, ip(y + 1) = tp(y). In this case, 
for <p(y) = 1 the integral I± describes a g-Meijer function, which was investigated 
by Slater in (S3- 

For |g| > 1, the equation A(j/+1) = R(q v )A(y) has the following general solution 

that is we have an effective q — > replacement and a reshuffling of parameters 
in lETI) . 

We remind the reader that q = e 27r4 °\ The parameter a gives a second scale, 
which may be used for generating a natural additional restriction upon A(y). The 
function q v is periodic under the shift y — > j/ + cr _1 , and (|2.7() satisfies the equation 
A(y + <7 _1 )/A(y) = x 1 / 17 ip{y + a^ 1 )/ip(y). We can fix by demanding that 

where R is another rational function of its argument. In accordance with the 
periodicity condition ip(y + 1) = tp(y), we have 



2niy\ _ 



R(e 2my ) 



where ij and tUj are arbitrary new parameters. Note that we cannot multiply the 
function R by terms like pe 2vlkv , k £ Z, p £ C, if they are different from 1, because 
then the periodicity condition for tp(y) will be broken. For \q\ < I, the general 
meromorphic solution of the difference equation for tp(y) is as follows: 

nJi(9'i7 1 e 2 -'»;9)~n;' =1 (tie- M »;'!W 

where g = er 2m l a is the modular partner of q. Indeed, for Im(<r) > we have 
Im(<r -1 ) < 0, and 1(2.8(1 is well defined. The function (p(y) in (|2.8() is an arbitrary 
elliptic function with periods 1 and <r _1 . It is characterized uniquely by the position 
of its poles and zeros in the fundamental parallelogram of periods containing 2k — 1 
free parameters, where k is the order of (p{y). Thus, the space of solutions is not too 
large: it becomes finite-dimensional (in the sense of the number of free parameters). 

Consider the regime |g| = 1. Denoting a = u>i/u>2 and assuming that Re(cr) > 0, 
we introduce the variable u = ywi- Now it is possible to choose the parameters 
tj,tj, etc in a special way, so that the infinite products (tjq v ;q)oo, (tje~ 2wty ; q)oo, 
etc in (12.7(1 and 1(2. 8|l combine into the double sine functions S(u + gf,w\,U2) for 
some c/j, where 

S(u;u 1 ,uj 2 )= )\ . , , (2.9) 

is a well defined function for \q\ — > 1. Indeed, it can be checked that the zeros 
and poles of 1(2. 9|) coincide with the zeros and poles of the function ^(wi + u>2 — 
u;uj)/T2(u;u>), which is a well-defined meromorphic function of u for oj\/u>2 > 0. 

In this case a is real, and if it is incommensurate with 1, then (p(y) = 1 (i.e., the 
function A(y) is determined quite uniquely). For a description of the properties of 
the double sine function and some of its applications, see |JMI IKLSl IN Ul lR"u2) . In 
particular, the integrals introduced by Jimbo and Miwa in |JM| as solutions of some 
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g-difference equations at |g| = 1 provided the first examples of g- hypergeometric 
integrals for q on the unit circle. Faddeev's concept of the modular double for 
quantum groups (see 0) is also related to the function (|2.9I) . 

Thus, the world of g-Meijer functions appears to be reacher than in the plain 
hypergeometric case. The introduction of the additional equation involving shifts 
by er -1 brought some new non-trivial structures in the integrals and reduced the 
functional freedom in the definition of meromorphic function A(y) to an elliptic 
function (p(y) containing a finite number of free parameters. 

Now, we turn to the single variable elliptic hypergeometric integrals. The general 
elliptic function of order r + 1 can be factorized as follows WW : 

J-J- 9 1 (« J +y;a,r) 9(wq<P, . . . ,w r qV;p) 



3=0 

9{t ,...,t k ;p) = l[8(t l ;p), 



k 



i=0 

where p = e 27 ™ T , Im(r) > 0, q = e 27rl<T . The parameter 7 is an arbitrary complex 
number, but i, = q Ui , Wi = q Vi satisfy the balancing constraint 

r r r 

^2{ Ui -Vi)=0, or = Il^i, (2.11) 

i=0 i=0 i=0 

which guarantees that the meromorphic function h(y) is doubly periodic: 

h(y + a' 1 ) = h(y), h(y + to' 1 ) = h(y). 

For t = a (which requires that Im(er) > 0), the function h(y) gives an explicit form 
of (p{y) in Q2SJ. 

In order to find the integrand A(y), it is necessary to solve the first order differ- 
ence equation 

A(y + l) = h(y)A(y) (2.12) 

in the class of meromorphic functions. The theory of such equations was developed 
long ago (see, e.g., |fja2| ). Obviously, since h(y) is factorized into the ratio of 
products of theta functions, it suffices to find a meromorphic solution of the equation 

f(y + l) = e(qV;p)f(y), (2.13) 

which leads to various elliptic gamma functions Q]. The simplest such function 
(|1.1U|I is defined from equation (|2.13|) only up to a periodic function <p(y + l) = (p{y) 
and, moreover, it requires that Im(cr) > (or |g| < 1), which was not assumed in 

We introduce the variable z = q y , so that the shift y — * y + 1 becomes equivalent 
to the multiplication z — > qz. Then the general solution of 12.12fl looks like this: 

where the balancing condition l|2.11|l is assumed and <p(y + 1) = tp(y) is an arbitrary 
periodic function. Using the reflection formulas 

T{pz,qz^ 1 ;q,p) = T{qz,pz^ 1 \q,p) 

= r(pg2,2 _1 ;g,p) = r(z,pgz _1 ; q,p) = 1, (2.15) 
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in l|2.14(l we can replace several elliptic gamma functions containing in the argu- 
ments z by those with arguments containing z . After that, A(y) would look 
closer to the integrands for the plain or g-Meijer functions, but in the elliptic case 
this does not increase generality because of the right-hand side of (|2.15() is trivial. 

In the region Im(cr) < 0, that is, for |g| > 1, the general solution of i|2.12|l can 
be written in the form 

Effectively, we have a permutation of parameters and a simple q — > q^ 1 substitution 
in the elliptic gamma functions in (|2.14|) (cf. the definition of this function for |g| > 1 
given in |FVj 1. 

Let us take <p(y) = 1. Then the function i|2.14[l satisfies two simple difference 
equations of the first order: 

A( 2/ + a- 1 ) = e ^A(y), (2.17) 
A(y + ra- 1 ) = f[ ( 2 - 18 ) 

Suppose that 1, er , tct -1 are pairwise incommensurate. Then the system of three 
equations l|2.12ll . I|2.17|l . and l|2.18|l determines A(y) uniquely up to a factor. As in 
the g-hypergeometric case, we can generalize equations (12.1711 and (|2.18(l . use them 
as natural tools for fixing the functional freedom in A(y), and get qualitatively 
different elliptic hypergeometric integrals in this way. 

The ratio A(y + ra^ 1 )/ A(y) in l|2.18(l is an elliptic function with periods 1 and 
(j -1 . Therefore, it is natural to demand that A(y + cr _1 )/A(y) be also an elliptic 
function with periods that, by symmetry, are equal to 1 and tct -1 . 

Theorem 2. Suppose that A(y) satisfies equation and that 1, rr -1 , rer -1 are 

pairwise incommensurate. Denote q — e - 2lTl /< 7 ; p = e 27rir/<r _ p or s i m pH c ity ; assume 
that Im(a) > (i.e., \q\ < I). If A(y + ra' 1 ) / A(y) is an elliptic function with 
periods 1 andrcr~ 1 , then for Im(r j 'a) > the most general form of the meromorphic 
function A(y) is as follows: 

where Y[ T j=a tj w j = Ef/Lo ^j^J ~ ^' ^ or I m { T /°~) <0, we have 

Proof. First, observe that for Im(cr) > we have Im(tr _1 ) < 0, automatically, that 
is \q\ < 1. Therefore, for Im(T/cr) > the function T(z; q,p) is well defined. 

The function A(y) in (|2.14ll gives the general solution of equation l|2.12[) . Suppose 
that A(y + a~ 1 )/A(y) is an elliptic function of order n + 1 with periods 1, to~ x . 
For Imfr/cr) > 0, this demand is equivalent to the following equation for tp(y): 

+ Q fr Q{t 3 e-^y-p) 
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where nj=o^J% = 1- Note that we cannot multiply the right-hand side of this 
equation by any constant different from 1, since this would violate the condition 
ip(y + 1) = ip(y). 

The meromorphic solution of (|2.21l) is 

p* 1 (Wje ^ ly ;q,p) 

where (p{y) is an elliptic function with periods 1 and a~ 1 . We can write 

-pV e(a e-^y-q) = ™ r(a J -e- 2 "^,6 3 e- 2 ^;g,p) 
7=i 6{bje- 2niy ;q) fJi r(a,je- 2 ™y ,bje- 2 ™vp; q,p) ' 

where Iljli a j^J = Therefore, we can absorb the function (p{y) into the ratio 
of elliptic gamma functions in (|2.22JI by changing n — > n + 2m and identifying 
4 = pa k ,w k = a k for k = n + 1, . . . , n + m and t k = b kl w k = pb k for k = 
n + m + 1, . . . , n + 2m. Since n, £j, tDj are arbitrary, without loss of generality we 
can set <p(y) = 1, which yields the desired expression Ij2.19|l . 
The function l|2.19|l satisfies the following equations: 



A( y + ^) = e^n ^, 2 ^ A(,), (2.23) 

A(y + r*->) = ft ff^ TT ^T^/ ? A(y). (2.24) 

The elliptic functions defined by the products JYj=o an< ^ IIj=o m W-'^ty have dif- 
ferent forms though both have periods 1 and er -1 . They are related to each other 
by the modular transformation a — > — I /a for the corresponding theta functions. 

Now, we consider the region Im(r/cr) < 0. Equations l|2.12|l and (|2.24|) are well 
defined in this case. They can be used for the determination of A(y), and it can be 
checked that, indeed, the function (|2.20() provides their general solution. Equation 
(|2.23fl is replaced now by the following one: 

A(») jiOfoe-^p-i)' [ b) 

that is, p in H2.21|) is changed to p^ 1 , and the parameters tj,Wj are replaced by 
p~ 1 ibj and p _1 tj) respectively. Using (|2.16() . it is easy to construct A(y) satisfying 
(|57T5| . J223), and in the \q\ > 1 region as well. □ 

In order to be able to work with q on the unit circle \q\ = 1, we need another 
elliptic gamma function: a kind of the elliptic analog of the double sine function 
Denote 

q = e "2 , q — e »i, 

p = e "2 ( p = e , (2.26) 

where cj, are some complex numbers. 

Suppose that uj\/w2 > and Im(o;3/w2) > (i.e., |p| < 1). Then we have 
Im(u!3 / u!i) = (uj2 / uJi)lm.(uj3 / U2) > 0, that is, \p\ < 1 automatically. Therefore, in 
the analysis of equation Ij2.12|l for \q\ = 1 it is necessary to assume that \p\, \p\ < 1. 
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Definition. Let \q\, \p\, \p\ < 1. Then we define a new elliptic gamma function by 
the formula 

In the limit p — > taken in such a way that, simultaneously, p — > 0, we get 
G(it; wi, o> 2 , w 3 ) - * 5' _1 (?i; u;i, w 2 ), where the double sine function S(u;uj) is fixed 
in (E2J). 

The function G(u;cj) satisfies the following three difference equations: 

G(u + wi;w) = % 27 ™^;p)G(u;w), (2.28) 
G(ii + ^ 2 ;w) =% 2 "^;p)G( M ;w), (2.29) 

r,/ , s 0{e 2m ^\q) rt{ 
G{u + uj 3 ; u) = G{u; oj) 

9{e ^q;q) 

— S(u; uji, u)2)S{<jj\ +0J2 - u; ojx, w 2 )G(u; uj). (2.30) 

For pairwise incommensurate u)\, w 2 , 0J3, these equations determine the meromor- 
phic function G(u;lo) uniquely up to multiplication by a constant, which follows 
from the nonexistence of triply periodic functions. 

The first equation requires that \p\ < 1, the second one requires \p\ < 1, and 
both of them do not impose any constraint upon q. The third equation (|2.30(1 
involves only the function S(u; uj\ 1 w 2 ), which is well defined for LO\ju)2 > 0, that 
is, \q\ = \q\ = 1. This means that the function G{u;uj) may be well defined in this 
unit circle region as well. 

In essence, the original elliptic gamma function Hl.lOfl has the same properties 
as the function (|1 . 1 ip . In a similar way, the function Ij2.27|l can be expressed as 
the following combination of the Barnes r 3 -functions up to an exponential of some 
polynomial in u of the third degree: 

T 3 (u; Uj, a> 2 , ^3)r 3 (^3 -ui-u; -ujj, -oj 2 , lo 3 ) 
r 3 (wi + oj 3 - u;wi,a>2,w 3 )r 3 (u - uji - uj 2 ; -0J1, -w 2 , coa) 

x r 2 (w 3 -oj 2 -u; -ll> 2 , uj 3 )T 2 (uj 3 - u;u>x,oj 3 ). (2-31) 

From this representation it follows that, indeed, (|2.27|) is well defined for real u>i, w 2 
with uji/lo2 > (and any complex w 3 ), like in the double sine function case. A more 
detailed analysis of this correspondence and an investigation of other properties 
of the function G(u; uj) will be given elsewhere. In particular, it is expected that 
G(u; uj) is the key function for an elliptic extension of the modular doubling principle 
for g-deformed algebras [FllKLS| . 

As a result, for \q\ = 1 we get a solution A(y) of equation (|2.12() by the mere 
replacement of T(q y ;q,p) in l|2.14() by G(yuji\uj). In the rest of this paper we limit 
ourselves to the case where \q\ < 1. Note that the region \p\ — 1 is not well defined 
in the elliptic functions setting. In a sense, the region of real w 3 /w 2 is reachable 
only at the level of the original Barnes multiple gamma functions. 

3. A THETA ANALOG OF THE MEIJER FUNCTION 

The integral corresponding to (|2.14(l may be considered as a kind of an elliptic 
extension of a particular Meijer function. The general Jacobi theta functions analog 
of the Meijer function appears in the case where h(y) is a quasiperiodic function 
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corresponding to the fourth case of the definition given at the beginning of the 
previous section, see 1)2.5(1 . 

Let Pz{y) — y^.-— i ony 1 be an arbitrary polynomial of the third degree obeying 
the property PzifS) = 0. The function defined by the integral 

cr r (',«,,,)-[ "^T^V ^ f 3 - 1 ' 

V w / Jc\[k=Q T VWkq v ;q,p) 

where C is some contour on the complex plane, may be called a theta analog of the 
Meijer function whenever the integral is well defined. Note that no constraints are 
imposed in (|3.1|) upon the integers r, s and the complex parameters tj , uik 
We have the following equation for the integrand A(y) of l|3.1|) : 

(3.2, 

&(w qy,...,w r q v ;p) 

where Pa(y) = -PaO/ + 1) — -^(j/) i s a polynomial in y of the second degree. From 
the considerations of [55] it follows that this h(y) is the most general function 
such that h is meromorphic in y (together with its inverse l/h(y)) and satisfies the 
quasiperiodicity conditions 

h(y + a' 1 ) = e ay+b h(y), h{y + tcj- 1 ) = e cy+d h(y) (3.3) 

for some constants a, 6, c, d. The function h(y) may also be interpreted as a general 
meromorphic modular Jacobi form in the sense of Eichler and Zagier |EZj . 

However, the integral (|3.1|l is not the most general integral leading to l|3.2|l . 
Using appropriate modifications of the integrands (|2.19() and (|2.20(l and replacing 
y by y /oj\, we arrive at the general theta analog of the Meijer function. 

Definition. In the definitions i|2.26[l of the bases, assume that \q\, \p\ < 1. Then, 
for \p\ < 1, the integral 



t,t 

a; uj 

w, w 



nLo Ffae g ; g ' P) ^ Ffee ' I; g "' g eftfa) dy (3.4) 

is called the general theta hypergeometric integral of one variable whenever it is 
well defined. For \p\ > 1, we set 

= /• in „r^— ^ii;.,"^ — , :g.P ■ v ,„^ (35) 

7c nLo ; g,p) Il" = o rfee-^-i;,-,?- 1 ) 

There are no constraints upon the integers r, s, n, m £ N and the complex parame- 
ters tj,tj,Wk,Wk- 

Both integrands of l|3.4|l and (|3.5I) satisfy the equations A(y + u>i) / A(y) = hi(y), 
2 = 1,2, 3, where hi are some quasiperiodic functions: hi(y + uJk) = e aikV+bik h(y), 
i 7^ k, with aik,bik being some constants related to the parameters t,t,w,w, a 
and uj. The integral l|3.5|l was determined by the condition that it has the same 
functions hi{y),hz(y) as (13. 4|) . For a special choice of parameters t, t, w, w, a, in 



THETA HYPERGEOMETRIC INTEGRALS 



13 



the limits \p\, \p\ — » or \p\ 1 — > the function G; s "j f w^' Q! ' a; ) ^ s reduced to 
the general g-hypergeometric integral considered in the previous section, see (|2.7|) 
and iJUJl. 

The general single variable theta hypergeometric series is defined by the following 
formula IS5I: 



/ to,...,t s 
s +iEr I ... .. ;a;q,p 



\w 1 , ...,w r 

6 (to, h, . . . , t s ;p; q) n 



DC 



uyL ,Ll,...,L s ,p,q )n e p 3(n) 



^ 9(q,Wi, . . . ,w r ;p; g) n 

Actually, these series are slightly more general than those introduced in |S5| . be- 
cause in that paper we considered only the case where a 3 = 0, but the generalization 
to (|3.6I) is straightforward. We note that the presence of cubics of the independent 
variable y in (|3.1|) or n in (|3.t)|) is natural since we are working at the level of the 
Barnes multiple gamma function i|1.2fl of the third order. 

Writing l|3.6|) in the form of the sum X)^o c ™ w ith c o — 1 > we easily see that 
Cn+i/cn — h(n), where h(n) is given by l|3.2|l with wq = q and y — n. This 
coincidence is not artificial. Consider the sequence of poles of the integrand in (|3.1() 
located at y = yo + n, n G N, for some yo- We denote by kc„, cq = 1, the residues 
of these poles. As y — > yo + we have A(y) — > nc n /(y ~yo — n) + 0(1). Now it is 
not difficult to see that 

lira ^ + i) = £?±I = Um %)=% + n). 
y— >2/o+" A(yJ c n y— >yo+n 

In particular, this means that the sums of the residues in the integral (|3.1|) that 
appear from appropriate deformations of the contour C, form the theta hypergeo- 
metric series l|3.6|l for some choices of the parameters. 

In accordance with the classification introduced in |S5| , the elliptic hypergeomet- 
ric series correspond, by definition, to h(n) equal to an elliptic function of n. Such 
series are called also the balanced theta hypergeometric series. They are defined 
by the following constraints imposed upon l|3.6|l : 

r r 

s = r, a 3 = a 2 = 0, I ' 1 1 " - ( 3 - 7 ) 

Similarly, the integral l|3.1(l will be called the elliptic (or balanced theta) hypergeo- 
metric integral if the conditions l|3.7|l are satisfied. Evidently, in this case h(y) in 
(13.21) becomes an elliptic function of y. 

When h(y) is an elliptic function of y and of all parameters Uj,Vj (we remind 
the reader that tj — q Uj , Wj — q Vj ) , we call (|3.6|l and (|3.1|) the totally elliptic 
hypergeometric series and integrals, respectively. As was shown in |S5| . in addition 
to the balancing requirement, such a property imposes the following constraints on 
the parameters 

tjWj = p = const, j = 0, . . . , r, (3-8) 

which are known as the well-poisedness conditions in the theory of basic hypergeo- 
metric series |GR| . The explicit form of the integrand A(y) for well-poised balanced 
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thcta hypergeometric integrals is 

where we have denoted z — q y and 7 = ol\. The parameter p is redundant, it can 
be eliminated by the rescalings ti — > p x ^ 2 ti, z — > p~ 1 t 2 z, but we keep it for further 
needs. Observe that without loss of generality one of the parameters in ()3.1[) can 
be set equal to one by a shift of y. 

Without the balancing condition, a theta hypergeometric series r +\E r is said to 
be well-poised if the constraints l|3.8[l are valid with wq = q, and very-well-poised if, 
in addition to l|3.8() . we have 

, _ ,1/2 , _ ,1/2 

tr-3 — to ^-2 ~~ — f 9' 

t r - 1 =t 1 /2 qp-W t r = -tl /2 qp^. (3.10) 
Such series take a simpler form 

t ,ti,...,t r -4,qt ,—qt ,qp 1 t ,-qp ' t 

,1/2 ,1/2 „ 1/2 ,l/2 „-l/2. 



_1 r + /+ * +1/2 +1/2 1/2 + 1/2 _ 1/2 ,l/2; a ; ( ?'-P 

\qto/t 1 ,...,qto/t r -i,t »-*o >P 7 V >~i> 7 V 



_ v 2 ^ e(t q 2n - lP ) r y\ 0{t m \P\q)n , n X„P 3 (n) 

^ e(t 05 p) ii %W^;p;g)n l 9j ' 1 j 

The essence of (|3.1U|I consists in the replacement of the product of four 9(tiz;p) by 
one theta function oc 8(t q 2 z 2 ;p) (this corresponds to doubling the argument of the 
#i-function). Very- well-poised series play a distinguished role in applications, in 
particular, they admit an appropriate generalization of the Bailey chains technique 
of generating infinite sequences of summation or transformation formulae S6 . 

In the case of integrals, we call (13.111 the very -well- poised theta hypergeometric 
integral if, in addition to conditions l|3.8|l . eight parameters ti are fixed in the 
following way: 

(* m _ 8l . . . , t m _r) = (±(pq) 1/2 , ±q 1/2 p, ±p 1/2 q, ±pq). (3.12) 
These constraints lead to squaring the argument of the elliptic gamma function 

m— 1 

TT F(tjz;q,p) = T(pqz 2 ;q,p) = — — (3.13) 

1 (z , q, p) 

(such a relation was used already in |S2| in the derivation of the elliptic beta in- 
tegral 11.1211 ). As a result, the integrand of the very- well-poised balanced theta 
hypergeometric integral takes the form 

AM-ff r " J ' ; "' ^"-^^'^y . (3.14) 

After imposing conditions (|3.12|) . the parameter p is no longer redundant, and its 
choice plays an important role. If we fix it as p — pq, then A(y) takes a particularly 
symmetric form 

A(y) = 'V y 2 \ \ ; \ e™, (3.15) 
l(z 2 ,z z ,Az,Az L ;q,p) 
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where A — (pq) Y[ jTo 9 tj ■ Clearly, the cases where m is odd or even differ from 
each other in a qualitative way. The choice m = 13 gives the simplest expression for 
A(y) and plays a distinguished role. Other simple choices, m — 9 or 11, correspond 
to particular subcases of the situation with m = 13. For m = 13 and 7 = we get 
the integrand of the elliptic beta integral 11.13|) . that is the simplest very- well- 
poised elliptic hypergeometric integral turns out to be exactly calculable when C 
is taken to be a special cycle corresponding to the unit circle on the z-plane. 

The sums of the residues of the function (|3.15|l for m = 13 are expected to form 
a 14-E13 theta hypergeometric series. However, the very-well-poisedness condition 
(I3.12|) results in the cancellation of theta functions in the corresponding ratios of 
the series coefficients h(n). As a result, we get only a lo-Eg very- well-poised elliptic 
hypergeometric series which, for 7 = 0, corresponds to the left-hand side of (| A. 11|> 
or (|6.8|l at n = 1 (for more details, see |DSlllS5] ). This shift of indices m — » m — 4 
brings in one more intriguing point related to the origins of the very-well-poisedness 
condition. It is necessary to find some deeper algebraic geometry explanations of 
the fact that in the single variable case "the nice things" (summation or integration 
formulae) are related to the number 14, the order of the initial elliptic function h(y). 
As was shown in |S5MS7| . in the multiple case this order raises to higher numbers, 
but in quite an intriguing way as well. 

Since T^pq) 1 / 2 z, (pq) 1 / 2 z -1 ; q,p) — 1, we may drop two parameters ±(pq) x / 2 in 
the condition l|3.12|l . For p = pq, 7 = 0, this yields 

TTl'-n 7 tjZ -1 : q, p) 

A (V) = J 2 -2 A A 1 V ( 3 ' 16 ) 

1 (z , z Z ,—Az,—Az ,q,p) 

where A = (pq) 2 Y\™L$ tj. For m = 11 this expression looks similar to i|1.13fl . 
but the different sign in front of A changes the things drastically, and it is not 
known whether the corresponding integral gets any closed form expression. 

In a more general setting, we can impose balancing and very-well-poisedness 
conditions upon general theta hypergeometric integrals (|3.4|) and 13.5fl . However, 
at the moment it is not known whether the simplest integrals appearing in this way 
admit exact evaluation. 

As far as the multivariable integrals of hypergeometric type are concerned, the 
general form of A(y) in the plain and q- hypergeometric cases can be deduced from 
the Ore-Sato theorem for Horn's series (see, e.g., |G(tR) for a detailed discussion). 
The general form of the multiple elliptic hypergeometric series or integrals is not 
established yet. We formulate it as an open problem — to find an elliptic or general 
theta functions analog of the Ore-Sato characterization theorem. In the following 
sections we give a series of examples of multivariable extensions of the very-well- 
poised balanced theta hypergeometric integrals associated with the root systems 
A n and C„. 

As to the further possible generalizations, it is natural to consider integrals of 
hypergeometric type for arbitrary algebraic curves or general Abelian varieties. 
Both would involve Riemann theta functions of many variables, appropriate gener- 
alizations of gamma functions and theta hypergeometric series. Some preliminary 
discussion of ideas in this direction can be found in jS7| . in particular, a special 
subcase of the s^7 Jackson summation formula was generalized there to Riemann 
surfaces of arbitrary genus. 
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4. Known C n elliptic beta integrals 

The following multivariable generalization of the Euler beta integral I|1.1J) has 
been introduced by Selberg | A AR| : 

[ ■■■ [ H xf-^l-Xjf- 1 Y[ \ Xj -x k \*i dxi---dx n 

Jo ^° l<j<n l<j<k<n 

tj r(a + (j - i) 7 )r(/3 + U - i) 7 )r(i + 31) , 4n 

liin r(a + /3+(n + i-2) 7 )r(l+ 7 ) ' l ' j 

where Re(a), Re(/3) > 0, and Re (7) > - min(l/n, Re(a)/(n - l),Re(/3)/(n - 1)). 
This integral has found many important applications in mathematical physics. 

At the first glance, the integrand of l|4. ljl does not fit the definition of the hyper- 
geometric integrals introduced in the preceding section. However, we can rescale 
y% — * Ui£ in <|2.1|) . and choose rational functions Ri(y) appropriately so that the 
limit e — > becomes well defined and yields a system of linear differential equations 
of the first order. As a result, we get dA(y) /dyi = i?j(y)A(y) and, apparently, 
(|4.1(l satisfies these conditions. 

Two types of multidimensional generalizations of the elliptic beta integral (|1.12|) 
to the root system C n were proposed by van Diejen and the author in |DS11 DS2 . 
One of them reduces in a special limit to the Selberg integral (|4.f () . Here we 
describe these elliptic Selberg integrals explicitly. For brevity, we drop the bases 
p, q in the notation for elliptic gamma functions from now on. We introduce the 
Type I integrand as 

A 7 (z;C*„) = —J— ]J T ~ 1 (zj*k,Zj z k 1 > z 7 lz l" z 7 lz k 1 ) 

^ ' l<j<k<n 

n -i-r2n+2 p/, , — 1\ 

TT jjf=g n T rZj,t rZj ) 

n} inz *, Z r*,A Zj ,A z ry ^ 

where t r G C, r = 0, . . . , 2n + 2, are free parameters and A = Ilr=o 2 The Type 
II integrand has the form 

jj 1 r (tZj z k ,tz j z- 1 ,tz7 1 ZkM7 1 z^ 1 ) 

^ j M"i</<L T{ ZjZk , ZjZ - k \ Z J^ Zi ^) 

n Wr=0 F(t r Zj, t r Z- ) 
x T{ z i z f,B Zj , Bz fy ( - 6) 



x 

3- 

where t,t r € C, r = 0, . . . ,4, are free parameters and B = t 2n ~ 2 Ylt—o By ^ wc 
denote the unit circle with positive orientation. 

Consider the first type of the C„ multivariable elliptic beta integral. We take 
\p\, \q\ < 1 and |t r | < 1, r = 0, . . . , 2n + 2, and assume that \pq\ < \A\. Then 



JT" z \ z n 

= 2 n nl no< r <,< 2 » + 2 r (^ 

(p;p)S, («;?)£, n^o'r^ 1 ^) 



(4.4) 
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The second type of the multiple elliptic beta integral has the following form. We 
take \p\, \q\, \t\ < 1 and \t r \ < 1, r = 0, . . . , 4, and assume that \pq\ < \B\. Then 

JT 71 z 1 z n 

(p;p)Scfeg)So/J[ r(t) rir=o n* 1 ^ x £) 

The elliptic Selberg integral of Type II 14.5(1 can be deduced from the Type I 
integral with the help of some interesting trick DS2 . The type I integral (|4.2() 
was proved in DS2] under a vanishing hypothesis, namely, that its left-hand side 
vanishes on the hypersurface of parameters A = ti n +i after an appropriate de- 
formation of T to an integration contour C that separates the sequences of poles 
in z converging to zero from those diverging to infinity. As p — > 0, both Type I 
and II integrals are reduced to Gustafson's well-known q-Selberg integrals |G11IG2| . 
which are related (for particular choices of parameters) to the Macdonald-Morris 
constant term identities and the Macdonald polynomials for various root systems 
including the Koornwinder polynomials 

As was shown in DST1 IDS3] . the sums of residues of the functions 1(4.2(1 and (|4.3|) 
form some multiple elliptic hypergeometric series. In particular, the multivariable 
10-E9 sum conjectured by Warnaar in |Waj can be deduced from 1)4.5(1 . The residue 
calculus for ((4.4(1 yields an elliptic extension of the basic hypergeometric series sum- 
mation formula proved in |DG| and [ML . Recursive proofs of these multivariable 
Frenkel-Turaev sums were given by Rosengren in [Roll IRo2] . The property of well- 
poisedness (or total ellipticity, in the case of elliptic hypergeometric series) plays an 
important role in such summation formulas. First examples of (plain) multiple hy- 
pergeometric series well-poised on classical groups were introduced by Biedenharn, 
Holman, and Louck in HBL . 



5. A NEW C n INTEGRATION FORMULA 

We define 

AlII(z;Cn) ^(2^ II ^( Zi zr\z^;p) 

V ' l<i<j<n 

11 11 y(z 2v z"A) ' [ ' 

i=lv=±l L[Z * ' 1 > 

where A = ttit 2 hq n ~ l . 

Theorem 3. Impose the following restrictions upon parameters: \xi\, \tk\ < 1, 
\t\ < |arj|, where i = 1, . . . ,n, k = 1, 2, 3, and \pq\ < \A\. Then 

f A /// / \ d>Z\ dz n 2 t r n t 1 / \ 

/ A (z;C„) = - — ; — II x 3 e{x l x j ,tx i x j ;p) 



VnU\T[ [ ^■l<r<s<3 T ( t r t sq' t ^ 1 f V (x t t k , tt k / X t ) \ 

{t} l\ { T(A/ Xl ,A Xl /t) 1 = 1 T(Aq^/t k ) ) ■ (5 ' 2) 



18 



V.P. SPIRIDONOV 



Proof. Consider the determinant 

det ( / ^{z^.Mq^' Mq 3 ~\h,txi 1 )— ) 

l<i,j<n \J T Z J 

= TT^f —■■■ — f[G l (z l )G l (z^)D(z), (5.3) 



where Ae is the integrand of the elliptic beta integral (|1.13|) with an appropriate 
choice of the parameters. The expression standing on the right-hand side of (|5.3I) 
appears after taking the integral signs outside of the determinant symbol, so that 
we get a multiple integral with 

f^i / \ ^v^i^ij ^i^2> ^«^3? ^ii/^i) 

l{Zl) = r(zf, Zl A) ' 

D(z) = det (9(ziti,z~ 1 ti;p;q) n - j 6(z t t2 1 z^ 1 t 2 ;p;q) 3 -i) . 

l<i,j<n 

The determinant D(z) can be rewritten as follows 
n"=i 0(si*2, Zi 1 t 2 - 1 p; q) n -i 



D(z) 



x det 



6(z i t 1 ,z^ 1 t 1 ;p; q) r 



l<i,j<n \6(q 2 - n /zit 2 ,q 2 n Zi/t 2 ;p;q) n -] 

In |Wa| . Warnaar computed the following elliptic generalization of a Krattenthalcr 
determinant IKrl: 



det 



/ 9(aXi,ac/Xi;p;q) n - 



i<j,i<« V bc/Xi;p; q) n -j 

= a(")g(S) J] X i 9{X i XJ 1 i cXZ 1 Xj\p) 

l<i<j<n 

n e(b/a,abcq 2n - 2l ;p; q)^ 
; i 6{bX u bc/Xi;p;q) n -i ' [ ' ' 

Using this identity for X,- = Zi, a = ti, b = q 2 ~ n /t 2 , and c = 1, we find 

£>(*) = (tif^ftVG) J] z^ziz^zrh^ip) 

\<i<j<n 

n 

x n^^ 2_n /*i*2>*i9" +2 " 2< A2;p;?)i-i- (5.5) 

i=l 

As a result, the determinant l|5.3ll yields an expression proportional to the left-hand 
side of the C n multiple integral in question (|5.2H . 
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Now we substitute the result of computation of the elliptic beta integral (|1.12(l 
into the determinant Q5.3JI . This yields 

2 n 

l.h.s. of O = — — — 

T(x l ti,t 1 t/x i ,x l t 2 ,t 2 t/x i , Xit 3 , t, t^q"-' 1 ,tit 3 q n ~ l ,t 2 t 3 q l ~ 1 ,t 3 t/xi) 



x 



n 



T(A/x i} A/ta, nA/t, Aq 1 ^ /t 2 ,Aq l ~ n /ti) 



x det (9(xiti,tit/xi;p;q)n-j9(xit2,t 2 t/xi;p;q)j-i) . (5.6) 

l<i,j<n 

By 1)5.4(1 . the determinant in the last line takes the form 

(ii^i)(SV(S) [] x j e{x i /x j ,t/x i x j ) P ) 

l<2<_7<n 

n 

x II % 2 ~7*i*2i, t iq n+2 - 2i /t 2 ;p; ?),_!. 
i=i 

Equating the resulting expression in (|5.6|) with the right-hand side of l|5.3[) . we get 
the desired integral evaluation (|5.2() . □ 

This is the first nontrivial multiple elliptic beta integral with a complete proof. 
It is not symmetric in p and q, unlike all other cases considered in the present 
paper. This fact suggests that there should exist yet another integral of similar 
nature that would be symmetric in p, q. 

The above method of computation of the taken C n integral represents a next 
step in the logical development of applications of determinant formulas to multiple 
basic hypergeometric series; see, e.g., the paper GK of Gustafson and Kratten- 
thaler, which was followed by Schlosser |Scll ISc2j and Warnaar |Wa| . Similar 
considerations for computing some multiple g-hypergeometric integrals were given 
by Tarasov and Varchenko in |TVj . 

6. An elliptic beta integral for the A n root system 

In this section we conjecture a multiple elliptic beta integral for the A n root 
system, which will be used in the next section for derivation of other nontrivial A n 
integrals. 

Conjecture. Let z,, i = 1, . . . , n, tfc, k = 1, . . . , n + 1, and fj,j — l,...,n + 2, 
be independent complex variables (n is an arbitrary positive integer). We denote 

A=ir^it k ,B=u^fi,^ 

i ngj (ig m^) n;=i 2 m^, 

A (z; An) - M « ^ r( ^) nS r(AB Zk ) ' (6 ' 1} 

where z x z 2 ■ ■ ■ z n+ i = 1. 

Suppose that the parameters tk, fj satisfy the constraints |tfc|j|/j| < 1, \pq\ < 
\AB\. Then the following integration formula is conjectured to hold true: 

/ A ; (z;^..^ = r i^ 

x —-, -t. . (6.2) 

nSrt^in^rt/- 1 ^) 
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For n = 1 this conjecture is reduced to the elliptic beta integral (|1.12(l . For 
arbitrary n and p — 0, we get a Gustafson integral proved in |G1| . Let us show 
that the two sides of l|tj.2fl satisfy one and the same difference equation. 

Theorem 4. Let I n (t,f) denote either side of Then this function satisfies 

the q-difference equation 

g'^n'lf^ * w f) =/„M> (o.3) 

and its partner obtained by the permutation of q and p. 

Proof. Denote the function H6.1|) by A 7 (z; t%, . . . , t n +\\ A„). It is not difficult to see 
that 

A J (z;t 1; . . .,qt r , . . .,t n+1 ;A n ) _ 'yr 9(t r zZ 1 ;p) 
A'(z; t n+1 ;A n ) ~ 11 9(ABz k ;p) ' 

so that equation (|6.3(l for the left-hand side of (|6.2J) is satisfied if the following theta 
functions identity is fulfilled: 

^ 9{Bt r -p) "+ 1 6{ABt 3 ;p) "+ 1 9(t r zZ 1 ;p) _ 

2_, 9 (A;p) 11 6(t r tf;p) 11 6(ABz k ;p) ' 1 ' > 

For ii = 1 this identity is equivalent to the well-known relation for products of four 
theta functions 

9{xw, x/w,yz,y/ z;p) — 9(xz, x/z, yw, y/w;p) 

= yw~ 1 9(xy,x/y,wz,w/z;p) (6-5) 

and equation coincides with that used in |S2| for proving the integral l|1.12|l . 

For n > 1, identity (|6.4|l can be established with the help of the Liouville theorem, 
much as in the arguments presented in DS2]. A simpler proof follows from the 
general theta functions identity given in |WW| . As was shown by Rosengren in 
|Ro2| . that identity can be rewritten as the following generalized partial fractions 
expansion of a ratio of theta functions: 

-pr 9{t/b k ;p) = 9{ta 1 ---a n /a r b 1 ---b n -p) U"=i^r/bj;p) 
■}J ■ 9(t/a k ;p) ^9(t/a r ,ai---a n /b 1 ---b n ;p)Y\ a 3=^i. a r/afp) 1 

where a±- ■ -a n ^ 6i • • • b n . Here we replace n by n + 1 and substitute a k = tZ , 
b k = zZ , and t — AB. As a result, we get an identity which is seen to coincide 
with l|b.4f) due to the relation z\ ■ ■ ■ z n+ \ = 1. 

In a similar way, for the right-hand side of l|tj.2fl we get 

I n {ti,...,qt r ,...,t n+ i,i) _ 9(A;p) yj 9{t r fj]p) 

J„(t,f) ~ 9{t r B-p) 11 9(ABf^ ;p y 

and in this case equation H6.3J) becomes equivalent to the identity 

n;^ 1 9{ABt 3 -p) ^ mi 1 ^(h t1 ^) ■ 
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If we substitute here f n +2 = B/fi ■ ■ ■ f n +i and divide both sides by 9(ABf n + 2 ]p), 
then we get Ij6.fi|) with n replaced by n + 1 and with a,j = tj 1 ,bj — fj ,t = AB. □ 

The equation derived above works in the space of parameters tk, whereas in 
|(xl| Gustafson used an equation in the variables fj for proving the p = case 
of the integral (|6.2(l . It is of interest that the latter equation does not admit a 
straightforward elliptic generalization, namely, the corresponding partial fraction 
expansion cannot be lifted to the theta functions level. 

Another argument in favor of the validity of formula (|6.2I) consists in the fact 
that, via the residue calculus, it generates a multivariable 10-E9 elliptic hypergeo- 
metric series sum for the A n root system, which was proved by Rosengren in |Ro2| 
and which was considered independently by the author in |S7j . 

Theorem 5. The residue calculus for the integral \b\ ty) yields the following sum- 
mation formula: 

j=l,...,n 

" 9(t l tj 1 q- N ^;p;q) Xt " 9{tj;p-q) w 
X .W 9(qt it f;p;q) Xi lUit^^-q)^ 

Q(b, c;p; q)\ x \ -A- 9(dt J ,et j ;p;q)x j 
9(q/d,q/e- 1 p;q)\ X \ fJi 9{t J q/b,t j qj 'c;p;q) Xj 

_ 8(q/bd, q/cd;p\q)\ N \ A 9{t J q,t J q/bc\p\ q) N] ^ 
0(q/d,q/bcd;p;q)i N \ j-J^9(tjq/b,tjq/c;p;q) N / 

where \X\ = Xi H h A„, \N\ = Ni H h N n , and bcde = q i+ \ N \ . For n = 1 this 

is the Frenkel- Turaev sum |FT| , and for p — it is reduced to the Milne 's multiple 
8$7 sum for the A n root system jJVlil] . 



x 



Proof. We scale tj for i = 1, . . . ,n from the region \ti\ < 1 to \ti\ > 1, and keep 
l^n+iU/jl < 1 together with the condition \pq\ < \AB\. During this procedure, 
some poles of the integrand A 7 (z; A n ) in l|6.2|l go out of the unit disk and, on the 
contrary, some of them cross over T entering inside. The outgoing poles are located 
at the following points: Zk = {tiq Xi ,i = 1, . . . , n} for each k = 1, . . . , n, and the 
number of such poles is determined by the conditions \tiq Xi \ > 1. The ingoing poles 
correspond to the points Z\ ■ ■ ■ z n = {t^ q~ Xi , i = 1, . . . , n}. 

We denote by C a deformed contour of integration such that none of the poles 
mentioned above crosses over C during the change of parameters. By analyticity, 
the value of the integral i|tj.2|) is not changing when C replaces T, that is the right- 
hand side of l|tj.2[l remains the same. If we start to deform the contour C back to 
T, we start to pick up residues from the poles by the Cauchy theorem. As a result, 
the following formula arises: 

J c n Li Zk U Jt3 Li Zk 

where R n = A I (z;A n ), and Rj(zi, . . . ,Zj) for j < n are sums of the residues of 
A 7 (z;t4„) corresponding to the poles crossing C. 
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We shall not derive explicit expressions for all coefficients Rj as it was done for 
the C n integrals in |DS1I IPS3] . For our purposes, it suffices to pick up only the 
residues that diverge in the limits fj — » q~ N Hj , Nj E N, for all j = l,...,n 
simultaneously. First, consider the residues appearing from the poles Zj = tjq Xj , 
where Xj are some integers such that \tjq Xi \ > 1. Straightforward computations 
yield 

n 

B$» (A) = Yl lim (1-tj q x * zj 1 ) A 1 (z; A n ) 

j=l z i^ t 39 Aj 

\i [%J q ) lH =1 T(t k q>*DqM,D-i<rMtt 1 q->*) 



*,i=i 



X 



nLi ( m^m^\- Xk )u;^mw k 



k=l 



ir k=1 r(ABt k q^)r(ABD-^q-M) 

(_l)A fc(? A fc (A fc + l)/2 

{q;q)oc(p;p)ocO(q;p; q)\ k 



™ (_ 1 )A fc9 A fc (A fc + l)/2 

11 (n-n\.Jr l -W\^ft(n-r l -n\, ' 



where £> = A/t n +i. The factors T(fjtjq Xj ) provide the required divergence in the 
limits fj -» q^ N 'tj 1 . We write i?^(A) = n n A(\;A n ), where 

n r -x fM -r f -r f jitlnupmtlnhD- 1 ) 



l<i<j<n 



\ i#fc 

X (g; «)« (p; p)« FlLi r(ASi fc ) r(^fl-i) ' 
and after a chain of simplifying calculations, A (A; A n ) takes the form 

A(X;A n ) = q^ J[ l '* ' P) 

l<i<j<n 'Pi 

" / 6(t 3 D;p;g) M e^Dg^+^-p) yr 6{fjt k ;p;q) Xk 



i=i y^f^D-.p; q) w 0(tjD;p) ^ e(qtjH k ;p; q) Xk 
6(cD,qD/AB;p;q) lx \ 6{dt kl et k ;p;q) Xk 



9(qD/d,qD/e;p;q)\ X \ 6(qt k /c,ABt k ;p;q) Xk 
where we have denoted c = t n+ \,d = /„ + i,e = f n +2- Now we substitute /j 



-JV, 



t ■ 1 (which assumes that = g \ N \cde) in this expression and introduce the 
parameter b = q 1+ \ N \ / 'cde. As a result, the function A(\;A n ) becomes equal to 
the summand on the left-hand side of 1)6. 8fl after the transformations tj — ► tj/D, 
b^>b/D,c^>c/D,d^>Dd,e^De. 

Now, we find the total number of residues of such type. There is permutational 
symmetry between the variables Zi, . . . , z n . Therefore, there are n! ways to satisfy 
the equalities z k = tjq Xj using each tj only once. The residues of the other outgoing 
poles located at z k = {Uq Xi , i = 1, . . . , n}, where at least one t% enters twice, do not 
diverge at fj — > q^ Nj tj 1 for some j. 
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We pass to the ingoing poles. It is not difficult to verify that the residues of 
A I (z; A n ) for the poles located at z n = tJ 1 q~ Xj / z\ ■ ■ ■ z n -i for some fixed j (or, 
equivalently, for z n+ i = tjq Xj ) are equal to the residues for the poles at z n = tjq Xj . 
Among the remaining poles in the variables Z\, . . . , z n —i, we must consider only the 
outgoing ones since only they may diverge as /& — » q~ Nh tu with k = 1, . . . , n, k ^ 
j. There are n ways to fix the variable z^ for which we shall consider ingoing poles, 
there are n ways to fix the parameter tj in the equation z n+ i = tjq x \ and there 
are in — 1)! appropriate outgoing poles with the required residue divergence. As a 
result, the contribution of these combined ingoing and outgoing poles is equal to 
n 2 (n — 1)!, and the total number of diverging residues (|6.1U|I is equal to (n + 1)1. 
Roughly speaking, the incoming poles imitate the (n + l)st independent contour 
of integration over z n +\, which enters symmetrically with z±, . . . , z n , and there are 
(n + 1)! ways to order these variables in the residue calculus. 

As has already been mentioned, 1)6. 9J1 is equal to the right-hand side of 1)6. 2)1 . 
Now we divide both these expressions by (n + 1)!k„ and take the limits as fj — » 
q~ Nj tJ 1 , j = l,...,n. Since n n — > oo in this limit, only the residues considered 
above survive in 1)6. 9)1 . and their sum is given by the elliptic Milne series 1)6.8)1 . As 
to the right-hand side, we get 

r.h.s. of (l6~21) 9(q/bd,q/cd;p;q)i N] 

hm 



f j ^q- N Hj 1 (n+l)\Kn 9(qD/d,q/Dbcd;p;q)\ N \ 

rr QjqtjDiqtj/Dbcip; q) Nj 
J-Ji 9(qtj/c,qtj/b;p;q) Nj 

which coincides with the right-hand side of 16.8)1 after the appropriate changes of 
parameters indicated above. The theorem is proved. □ 

Formula ()6.2)) generates the following symmetry transformation for integrals: 

fr nsf- 1 ) r ngj ng ^ uu dz 3 / Zj 
Vr nssj 1 ) r mil ng n^\M uu ^hi f6n) 

Here t, fj, Sj,j — 1, . . . ,n + 2, are free independent variables, B = n^=i fj> & — 
n"=i Sj, and it is assumed that \t\, \fA, \ Sj \ < 1, \pq\ < \t n+1 B\, \t n+1 S\. In order 
to derive this identity, it is necessary to consider the 2n-tuple integral 

i r mil mil rite, ^i 1 ) KXU ntz^wA 



(2w)» A 2 „ 11" T{z iZ j l ,w iW J l ) Y\li\ T(t^Bz k , t^Sw- 1 ) 

dz\ dz n dwi dw n 
Zi z n wi w n 

where z\ ■ ■ ■ z n+ \ — w\ ■ ■ ■ w n +\ = 1. Integration with respect to the variables 
Zk with the help of 1)6.2)1 makes this expression proportional to the left-hand side 
of 1)6.11)1 (after the replacements Wk — > z^ 1 ). Changing the order of integration, 
which is allowed because the integrand is bounded on T, we arrive at the integral 
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standing on the right-hand side of l|6.11|l . up to some coefficient. After cancelling 
common factors, we get the required identity. For p — this reduces to the Dcnis- 
Gustafson transformation formula |DG| . which describes a Bailey transformation 
for a terminating multivariable lo^g series. It is natural to expect that in our case 
formula (|6.11() yields a Bailey transformation for some terminating A n multiple 
12-E11 series appearing from sums of residues of the corresponding integrals. 



7. Some other A n integrals 



Now from conjecture l|6.2|l we derive a number of different multiple ^4„-integrals. 
Denote 

M l<*i<n+l T &*j >*i Z i) 

'it 1 r (i i z 3 , t 2 Zj , t 3 Zj , U zj 1 ,t 5 zr 1 ) 
x =i — . (7.1) 



Theorem 6. Suppose the validity of the conjectured A n and C n multiple elliptic 
beta integrals if 6'. 2\) and \4-4\j > respectively. Then, the following two integration 
formulae are true. For odd n = 2m — 1, we have 



I A"(b;A,) 



dz\ dz n [n + 1)! 



zi z n (q;q)^(p;p) r c 

<3t 



r(t m ,s m , s^-Hits) ni<,<,< 3 r(t m - V,-) 



nL4 r (* 2m - 2 s m - 1 *i*2t 3 i fc ) 
x t^- 1 ms) j , vst-Htts) Ui< t<e <3 nv-^uu) 



For even n = 2m, we have 

A^(z;4„)^i...— 



U . . (izi <iz„ (n + 1)! 



Zl z n {q;q)^(p;p)^ 



TT 1 1 »'*)- ' I':. I I 

X 



5 

i— 

r((t s y,v s j-H A t 5 ) n- = i n*=4 m^y- 1 ^) 



^ nL 4 r^'-^^-i^ 1 n- =1 u) 



3=1 



Proof. The proofs follow the procedure used by Gustafson in |G2| for proving the 
p = cases of the integrals (|7.2ll and (|7.3I) . We start with the case of odd n = 2m— 1. 
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Consider the following (4m — l)-tuple integral: 

T 4m-1 n-7=i ; ^ r^z- 1 ) ru ±1 rii<i<i< m r «<> <™7vr^, zjvj 



TT T(zi(ts) m 2 IlLi f fc) TT dzj^ t~t / d/Wj^ dxj^ 

/ / nLir(<- 1/2 t fc ^)r(^r- 2 s - 1 /2t li2t 3 )n ^ 4r(s -i/2^ a; , ) ^ 



where YYi=i z j = 1- Using the exact C„ integration formula of type I (see (14.411 s ). 
first we take integrals in (|7.4I) with respect to the variables Wj,j = 1, . . . , to, and 
after that with respect to Xj, j — l,...,m. The resulting integral is equal to the 
left-hand side of (|7.2(l up to the factor 

2 2m (m!) 2 T(s-H 4 t 5 ) 

!<fi<3 r^ 1 ^) r(*»-v»- V nJL **) ' 

In this two step procedure, we need the following restrictions upon the parame- 
ters: 

|*| <1, |*1,2 >3 | < |f| V2 , \pq\ < \t m - 3/2 ht 2 t 3 \ 

and 

5 

\s\ < 1, |t 4)5 | < k| 1/2 , Iwl < \t m - 2 s m -V 2 l[t k \, 

k—l 

respectively. However, the resulting expression can be extended analytically to 
the region |t/-| < 1, k = 1,...,5, \pq\ < \{ts) 2m ~ 2 n*=i *fc| without changing the 
integral value. 

Since the integrand function in (|7.4(l is bounded on the unit circle, we can change 
the order of integrations. First, we take the integrals over Zj, i = 1, . . . , 2m— 1, using 
the A n -formula (|6.2[) . Then we apply formula l|4. 4f> in order to take the integrals 
over Xj , j — 1, ... , m. Finally, we apply the intrinsic elliptic Selberg integral (|4.5() 
for taking the integrals over Wj, j = 1, . . . , m; this leads to the following expression: 

{2m) 4m - 1 (2m)\(2 m m\) 2 T{t m , s m ^ -A- rjt^f^V QLi *«) 

(fertooCg;?)^) 4 " 1 - 1 r((fe) m ,t™s m - 1 i 4 i5) ^ r^-V"- 1 ^ 1 nf =x *0 

™ r((f S y , fy- 1 ^) ni<,< fc < 3 r(^v^*i*fc) nL 4 nli m a y-H k u) 
M nLi r((i S )«+i-2t- 1 n - =1 *o nL 4 r^-s^-v n s =i ^ 



X 



As a result, we get the needed integration formula for odd n — 2m — 1. 
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In order to prove l|7.3|l . we consider the 4?n-tuple integral 

t*- n-"=i; & n^j 1 ) ru ±1 n^^™ r««#, <^) 

x fi r(^(t s )^-inLi^) Ml™' ^ 



X 



v i| 1 \r(t m - 1 /2t 1 t 2 ^,a"-V2t 4 t 5 a;y ) «;2'' ,xf ) J ] W Z fc ' 



where z j = 1- Repeating the same trick as in the case of odd n (that is, 

integrating successively with respect to the variables Wj and Xj and then changing 
the order of integrations in this expression) , we get (|7.3|l . □ 



In a similar way, we can establish elliptic analogs of the A n basic hypergeometric 
integrals of Gustafson and Rakha GuR . 

Theorem 7. Suppose the validity of the A n and C n multiple elliptic beta integrals 
H\) and \4-4\I j respectively. Denote 



(7.6) 



where A — t n+2 JXi=i b an d j=i z j = 1 ■ Then the following two integration 
formulae are true. For odd n = 21 — 1, we have 



_ ni U-A\— dZn = ( n + i y- r ( f ''IIfc=i f fc) , 7 
T- n) Zl '-- z n fag)" fop)" T^rifcLiife) 

ntin =21+1 T{ttitj) ril<i<j<2i l<i<j<2i+3 

n-ii r(^+ v n'L + i 3 **) nS +1 r(*+v nS 3 **) 

For even n — 21, we have 

d Zl dz n (n + 1)! rdlS **, t* 4 " 8 **) 



Jm l ' " J zi'"z„ (g;g)^(p;p)^ r(i<+ 2 h 

Yii=i rij=2/+2 r(ttjtj) Ili<i<j<2/+i r(tt,tj) IIi=2;+2 r(t' +1 ti) 

n "il 1 r(^+2i"l n «+4 tfc) n «+4 +2 r(t Z + l t . rjW tk) 
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Proof. In accordance with the procedure used in |GuR| . we consider the (31 — 1)- 
tuple integral 

uZ n^zj 1 ) n % r(f n-io W) 

1 nS^tc^o 



X 



X 



u=±ll<i<j<l j=l l \ w 3 > l llfc=22 + l Z k w j I 

dw\ dwi dz\ dzii-\ 



Wl Wi Z\ Z 2 l-1 

where YlfLi z i = ^ an d = Ilfc^+i l k- Integrating with respect to the vari- 
ables Wj with the help of formula 14.4J) . we get the left-hand side of (|7.7[l up to 
some factor. Changing the order of integration, we can integrate over Zi using (IG.2I) 
(where it is necessary to change Zk to z^ 1 ) and then over Wj using 14. 4|) . Equating 
two expressions, we arrive at formula (|7.7|l . 

In a similar way, in the case of even n — 21 we consider the (3Z + l)-tuple integral 

nS 1 (nS T{t^ ZiWj ^ Zi wf) ns 1 r(^7 x ) 



l+l T-r2;+5 p/j.l/2-i 

w TT TT T-i-1/ v v v -v\TT 1 lk=2l+l L V L L-kWj) 

x 11 11 r (WiW^WiWj )[[— — — — 



X 



dzui du>i + i dz\ dz 2 i 
wi ' ' Wl+i Z\ ' z% ' 



where Hi^Ji 1 z « = 1 an d £20+5 = t l Ylk=i z k- Repeating the same trick as in the 
preceding case, we get the desired integration formula 1)7. 8|l . □ 

Sums of residues for the derived integrals H7.2fl - H7.8D f° rm elliptic hypergeometric 
series on the A n root system that differ from the series <|6.8ll introduced in |Ro21IS7j . 
We skip their consideration and formulate only a conjecture concerning the elliptic 
extension of Theorem 1.2 in |GuR| . 

Conjecture. Suppose that N is a positive integer and IlsUi z k — Q~ N ■ Then 



E 



nr J =i;* # ^(*i*7 l )A 4 -A 1 n?=if(* n+i *7 1 n2i 1 3 **)-A J 

Ai + ... + A„=Af 



rt is even, 



0(1)-N . ,, (7J) 

~i is odd, 



n£» 3 +1 e(t(»+D/^ l )-A.0(t("+ 3 )/ 2 nr= n 3 +i to-* ' 

where 6»(a) A = 6(a;p;q) X - 

These summation formulas are expected to follow from the residue calculus for 
the integrals 17.7J1 and (|7.8|) . Some evidence in favor of conjecture l|7.9|l is provided 
by the following theorem. 
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Theorem 8. Denote t = q 9 ,ti = q 9i ,i — 1, . . . ,n + 3 (so that Y^j=x 9j ' + N = 
0). The series c(A) standing on the left-hand side of |7.ff| ) is a totally elliptic 
hypergeometric series, that is, the ratios of successive series coefficients 

u <w c (^ 1 ' •• • , Afe + 1, . . . , A„) 

= 77 \ \ 

c(Ai, . . . , \ n ) 

are elliptic functions of all unconstrained variables in the set (Ai, . . . , A n , g,gx, ... , 
g n +3)- Moreover, the functions ft.fe(A) are SX(2,Z) modular invariant. The ratios 
of the expressions standing on the two sides of \7.iJ^ are elliptic functions of g 
and n + 2 free parameters in the set (<?i, . . . ,17,1+3), and these ratios are modular 
invariant as well. 

We skip the proof of this theorem, because it consists of quite long but straight- 
forward computations whose structure was described in detail in |S5I IS7| in the 
process of similar considerations for different elliptic hypergeometric scries summa- 
tion formulas. Using the fact that there are no cusp forms of weights below 12, as 
in |DS11 IS7| from this theorem we deduce that relations l|7.9|) are valid in the small 
a expansion up to the terms of order of a 12 . This gives also yet another example in 
favor of the general conjecture of |S5| that all totally elliptic hypergeometric series 
are automatically modular invariant. 

The integrals l|7.2|l and i|7.3[l are expected to generate A n summation formulas 
similar to (|7.9|l . It is also natural to expect that all multiple elliptic beta integrals 
described above lead to integral representations for various multiple \2E\1 elliptic 
hypergeometric series generalizing the single variable formula announced in |S4| 
(see Appendix B for the proof of it). We suppose that, sas in the A„-case, there 
exist several types of elliptic beta integrals and elliptic hypergeometric series sums 
associated with the D n root system (see, e.g., |Ro2l IS7| h but their consideration 
lies beyond the scope of the present paper. 

8. Relations to the generalized eigenvalue problems 

Consider the very-well-poised theta hypergeometric series (I3.11f) with the addi- 
tional constraint (— q) n e P3 ^ = (qx) n for some x G C. Special notation for such 
series was introduced in |S6| : 

r+iK(^o; ti, ■ ■ ■ , t r -4; QiP', x) 

^ 8{toq 2n ;p) tt 0(t m ;p;q) n 
For l|8.1[) . the balancing condition 1)3.7(1 is reduced to 

m— 1 

As p — > 0, the r +iV r series are reduced to r _iW r _2 very-well-poised g-hypergeome- 
tric series of the argument qx (in the notation of |GR| ) . 

The balanced 12 Vn series with x = 1 plays an important role in applications. 
For instance, the elliptic solutions of the Yang-Baxter equation derived by Date et 
al in jD-()ll ID-()2j are expressed in terms of such series for a particular choice of 
parameters |FT| . In all cases of the 12U11 function to be considered below, we have 
x = 1; therefore, we omit dependence on this unit argument from now on. 
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We denote £(t) = i2Vn(t ;ti, . . . ,tr;q,p), where ]lL=i ^™ = *o<? 2 , an d assume 
that this series terminates due to the condition t m = q~ n , n € N, for some m. In 
|SZ1I ISZ2] . the following two contiguous relations for £ (t) were derived: 

£(t)-£(t ;ii,...,t 5 ,<r 1 t 6 ,<#7) (8.2) 
O(qt ,q%,qt 7 /t 6 ,t 6 t r /qt ;p) A #(t r ;p) c/ 2 . . . . , s 



e(qt /t e ,q 2 t a /t 6 ,to/t 7 ,t 7 /qt ;p) ^ 6(qt /t r ;p) 

]^0(f r W#o;z>)£(9 2 *o;s*i, ■■■ ,qh,te,qh 



O(t 6 /qt ,t 6 /q 2 t ,t 6 /t 7 ;p) ^ 



]^6i(M7/(^o;p)£(g 2 io;gii, ■ ■ ■ ,qh,t 7 ) 



6(t 7 /qt ,t 7 /q 2 t (h t 7 /t 6 ;p) ^ 

5 



= x I[fl(g*0/*r;P)g(t)- (8-3) 

9{qt a ,q 2 t ;p) £± 

Their combination yields 
fl(tr ^/t r ,qto/tr,p) 0{qto/t6tr . p) pfot tm u q - %i qh) _ m) 

+ e(y fa,qto/t 6 ; P ) d( ^ A7tr ) (f (t t . . t ^ - lt7) _ f(t)) 

6{qU/t 7 ,t & /t 7 -p) 

5 

+ 9(qt /t 6 t 7 ;p) J] 0(£ r ;p) £(t) = 0. (8.4) 

r=l 

As p — > 0, these three equalities are reduced to the contiguous relations for the ter- 
minating very- well-poised balanced lo'&g series of Gupta and Masson |GM| . Similar 
contiguous relations at the level of a<& 7 functions were constructed earlier by Ismail 
and Rahman |IR| . 

We change the parametrization of the ^-function and consider relation 1)8. 4J for 
the following function: 

R„(z;q,p) = 12 Vn ,t 3 z,—,q , — ;q,p ), (8.5) 

\£4 £2^4 2 £4 / 

where ^4 = nr=o^»" We re pl ace the parameter to in (|8.4|l by £3/^4; the variables 
ti,t2, and t3 are replaced by q/tot4, q/tit^, and q/titi) the variables £4 and £5 
by and ^4g"~ 1 /t4; and the variables t§ and t 7 by £32 and £3/2, respectively. 
As a result, we see that R n (z;q,p) provides a particular solution of the following 
finite-difference equation: 

2V(*)=0, P m = ^(z)(T-1) + ^(z- 1 )(T- 1 -1) + Km , (8.6) 

where T is the g-shift operator, Tf(z) — f(qz), and 

? z 9 / 0(z z ,qz 2 -,p) 

3 
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The functions f(z) — R n (z; q,p) solve Ij8.fil) for fj, = q n ,n € N. 

Equation (j8.6j) looks like a nonstandard eigenvalue problem with the "spectral 
parameter" (i; indeed, it can be rewritten as the generalized eigenvalue problem 

V v f(z) = XV s f(z), (8.9) 

where the spectral parameter A is 

a= ) qt ;; v,p { (8.10) 

and the operators 2?f, are obtained from T>^ after the replacements of fi by 
arbitrary gauge parameters £,?y e C, £ ^ rjp k ,qt4p k /Arj, k £ Z. Application of the 
theta function identity ljfi.5j) to equation 1|8.9J) yields 

which shows that the gauge parameters £, ?y drop out completely from the equation 
determining f(z), T>^f(z) = 0. 

A three term recurrence relation for the functions R n (z;q,p) was derived in 
|SZ1I |SZ2 . It appears from formula (|8.4(l if there we replace t§ by q~ n and tj 
by Aq" -1 /^, and substitute t\ — > qjt§t^t2 — » <ljt\t±,tz — > q/t^ti, £4 — > ^z, 
^5 ~~ * ^s/- 2 - After some work, this foirmula can be represented in the form 

(7(2) - a n+ i)B(A9 n_1 /*4) (R n +i(z;q,p) - R n (z;q,p)) 
+ (7(2) - j3 n -x)B{q- n ) (R n ^(z;q,p) - R n (z;q,p)) 

+ 5 (7(2) - 7 (t 3 )) i?„(z; g,p) = 0, (8.11) 

where 



# It -*s_ ak. 9 g <? 2 ^ g 2 a . n 1 

w V t 4 x' t 4 x' t ti ' tot 2 ' M2 ' A ' A?y ' "J 

S W = / , 2 2, 2 \ ' ( 812 ) 



5 = 61 -a->r-r> r-r.r-r >*3»7, — ;p , (8.13) 
\ A t i4 iim t2M V J 

{x) = p^ML (8.14) 

a n = j(q n /U), (3 n = jiq^A). (8.15) 

Here £ and r\ 7^ £,P k ,£, 1 p k , k G Z, are arbitrary gauge parameters (they are not 
related to in the difference equation, but we use the same notation). Substi- 
tuting H8.12M8.14jl m l|8.11j) and applying identity Ij6.5j) . we see that the auxiliary 
gauge parameters £, r\ drop out completely from the resulting recurrence relation. 

Since B(q~ n ) = for n — 0, the indeterminate is not involved in (|8.11|) 
for n = 0. We can say that R n {z; q,p) are generated by the three term recurrence 
relation (|8.11|) for the initial conditions J?_i = 0, i?o = 1. All recurrence coeffi- 
cients in IjS.lljl depend linearly on the variable 7(2), which absorbs z-dependence. 
Therefore, R n {z; q,p) are rational functions of 7(2) with n being the degree of poly- 
nomials in 7(2) in the numerator and denominator of R n . Moreover, the poles of 
these functions are located at 7(2) = at, ■ ■ ■ , ot n . 
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For a particular choice of one of the parameters and a discretization of the 
values of z, the functions R n (z; q,p) yield elliptic generalizations of Wilson's finite- 
dimensional gF$ and lo'&g rational functions |Wij . They were derived in |SZ1| 
from the theory of self-similar solutions of nonlinear integrable discrete time chains 
(for a brief review of the corresponding approach to special functions, see |S1I 
IS3|). Discrete analogs of equations l|8. 6|) . (|8.9(l valid for the latter finite-dimensional 
system of functions were derived in [J3Z3 with the help of self-duality. An equation 
satisfied by io$9 functions, appearing from R n (z;q,p) in the p — > limit, was 
investigated by Rahman and Suslov in |RS1 . The general three term recurrence 
relations l|8.11[) were considered in (Zhj and, in a different form related to Ru 
continued fractions, in jlJVLj . 

The solutions of the generalized eigenvalue problems are known to be biorthogo- 
nal to each other; see, e.g., |SZlllSZ4llZh| and the references therein. Here we would 
like to demonstrate that the elliptic beta integral (|1.12|) serves as the biorthogonal- 
ity measure for solutions of equation l|8.9|) . Consider the scalar product 

f dz 

/ A E (z;t)^(z)(V ri -XV i )^z) — 1 (8.16) 
Jc z 

where A E {z\ t) is the integrand of i|1.12|l and Q(z), ^(z) are some complex functions. 
The expression (|8.16|) can be rewritten as 

dz 

A E (z;t){K v - V v (z) - y r) (z- 1 ))4'(^)$(z) — 

c z 

r d7 

+ / Ajsiq-^Vniq-^Wq-^Wz) — 
JC- z 

r r\ 7 

+ A £ ( g z;t)^(( Z - 1 z- 1 )fM$(z)--A{ ?? ^e}, (8.17) 
Jc + z 

where {rj — > £} means the preceding expression with rj replaced by £. The inte- 
gration contours C± are obtained from C after the scaling transformations z — > 
q ±1 z. Suppose that the poles of A^(z; t) and the singularities of the functions 
^(z), 1 i l (q ±1 z) do not lie in the region swept by the contours C± during their de- 
formations to C . Then 18.17|) takes the form 



/ A E (z;tMz) (2J£-A2>f )*(*)-, 
Jc z 

where the adjoint (or transposed) operator T>J has the form 



(8.18) 



,-i . 



T A E (qz;t) _ x _ x A E (q 1 z;t) _ x , 
^=A^ty^ 9 Z )T+ A E (z;t) V & Z)T 

-Vt(z)-Vs(z- l )+Kt. (8.19) 

Suppose that $>\(z) is a solution of the equation (T> v — XT>^) &(z) = and $>\'(z) 
solves the conjugate equation (v^ — A'"D^J ^( z ) = for some A'. Both these 
functions can be multiplied by arbitrary functions f(z) satisfying the condition of 
periodicity on the logarithmic scale, f(qz) = f(z). After the replacement of $>(z) 
and *ff(z) in (|8.16l) and (|8.18|) by &\(z) and ^\i(z), these expressions become equal 
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to zero. In particular, 18.18JI yields the relation 

Jc z 

= (A' - A) / A B («;t)* A («)lf*v(«)- = 0, (8.20) 
Jc z 

which shows that for A' ^ A the function $\{z) is orthogonal to Vj'fy \>(z). 
We find a function g{z) such that 

g- l (z) (Vl - XDf) g(z) = V v - XV £ . 

After substitution of the known expressions for Ae{qz] t)/Ag(x; t) and the spectral 
parameter (see A in (|8.10() '). we get the following equation for g{z): 

e(jL. ,2M Az,^;p] 

\ t 4 z' * 4 q 2 z>rl 

9{Qz) = — 7"i — r»(4 

u \^ t 4 ' *4Z ' qz ' <? 

which is solved easily: 

p/gp M_ Az £ ) 
r \ > t 4 ' tgz ' z ' in 22") 

V i 4 ' t 4 z' <j ' qz >^>r/ 

Here we have neglected the arbitrary factor /(gz) = /(z), which has already been 
mentioned. As a result, we get a direct relation between $a(z) and ^x(z): ^xiz) — 
g(z)Q\(z), where g(z) depends on A as well. 
We denote A„ = A| A1=g n and 

JnW = £/(z)U=«" 



Substituting <I>A n (;z) = R n (z]q,p) in the derived biorthogonality relations, we see 
that the functions R n (z;q,p) are formally orthogonal to T>J g m (z)R m (z; q,p) for 



The general considerations of SZ4, Zh show that Rn{z; q,p), which are rational 
functions of j(z) with the poles at j(z) = a%, . . . ,a n , are orthogonal to other 
rational functions of 7(2), which we denote as T m (z] q,p), with the poles at 7(2) = 
(3i,...,{3 n . The choice of a n ,f3 n and the other recurrence coefficients in (|8.11() 
determine R n and T„ uniquely, so that permutation of all a n with (3 n permutes R n 
and T n . In our case, we see that parameters f3 n are obtained from a n = r y(q n /t4) 
after the replacement of £4 by pq/A. Equivalently, this replacement converts j3 n to 
a n . An important point is that the weight function A^(z,t) is invariant under such 
a transformation. Therefore, we can get T„ out of R n simply by the £4 — > pq/A 
involution, which yields 

rpf \ rr f At 3 AAA t 3 _ n Ag"" 1 \ 

T n {z;q,p) = 12V11 ;—,—,—,t 3 z,—,q , — , (8.22) 

\ q t ti t 2 z U J 

where the dependence on p in the parameters drops out due to the total ellipticity 
of this 12 Vn series. Comparing with the previous consideration, we see that 

^j9n(z)R n (z;q,p) = p n T n (z;q,p) 

for some proportionality constants p n , which are of no importance for us here. 
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Thus, the operator formalism developed above leads to the following formal 
biorthogonality relation: 



for some constants h n . Suppose that C = T and \t r \ < 1, \qp\ < \A\. Some poles of 
the functions R n ,T m cancel with zeros of A^zjt). The remaining ones approach 
T as the indices n, m increase. Starting with sufficiently high values of n and m, 
the contour T stops to satisfy the conditions used in the derivation of l|8.18l) and 
must be deformed. As is shown in Appendix A, l|8.23[) is true if C separates the 
points z — to,i.2,3P a 9 h 7 t i p a q b ^ m 1 and A~ 1 p a+1 q b+1 ~ n , o,ieN, from their partners 
with the inverse z — > z^ 1 coordinates. 

Relation (jO^ seems to remain true even if we multiply R n {z; q,p) or T n (z; q,p) 
by an arbitrary function f{z) with the property f(qz) — f(z). However, such 
nontrivial f(z) must contain singularities which are crossed over when the contours 
C± are deformed to C (otherwise f(z) = const). Therefore, the influence of such 
additional factors should be considered more carefully. Moreover, only for very 
special f(z) the normalization constants h n may admit exact evaluation. 

The weight function Ae(z, t) is symmetric in q and p, whereas neither R n (z; q,p) 
nor T n (z; q,p) possess such a property. We can try to restore this symmetry using 
the freedom in the factor f(z) = f(qz). Take f(z) = Rf.(z;p,q), k £ N, that is, to 
the functions R n (z\q,p) themselves with the permuted bases q and p. Then, the 
product 



satisfies two generalized eigenvalue problems: (18.611 and the p-difference equation 
obtained from it by the permutation of q and p. For (|8.6|) we should have u = q n , 
and for its partner u = p k . The function (|8.1U|) does not change under the sub- 
stitution fi — > pn. Therefore, the choice fi = q n p k , n, k £ N, gives "the spectrum" 
for both generalized eigenvalue problems. The first factor of R n k(z) is a rational 
function of j(z;p) (we indicate the dependence on the base p explicitly), but the 
second is rational in 7(2; q). Therefore, for generic q,p it is necessary to view the 
functions R n f c (z) not as rational functions of some variable but as meromorphic 
functions of z. 

In the same way, the series termination condition tg = q~ n in (|8.5(l may be 
replaced by = q~ n p 1 , which terminates simultaneously the 12V11 series for 
Rk(z;p,q). The property of the total ellipticity of the balanced r +iV r (to', t\,..., 
t r -4,; q,p) series plays a crucial role at this place: any parameter t\, . . . , t r _5 may 
be multiplied by an arbitrary integer power of p without changes (note that the 
parameter to plays a distinguished role and the series are invariant under the trans- 
formation to —* P 2 to)- 

As a result of the doubling of eigenvalue problems, the functions R n j~(z) turn out 
to satisfy a quite unusual biorthogonality relation (|A.13|) characteristic of functions 
of two independent variables, which was announced in |S4j . A rigorous consideration 
of all these biorthogonalities associated with the function R n (z; q,p) with complete 
proofs is given in Appendix A. 

In fact, there is a deeper relationship between the structure of the elliptic beta 
integral (|1 . 12|) and the biorthogonal functions R n {z;q 7 p),T n (z;q,p) and the solu- 
tions of equation H8.6J) than it is indicated in this section. We hope to address this 




(8.23) 



Rnk(z) = Rn{z; q 7 p)Rk(z;p 7 q) 
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later on. As far as the multivariablc generalizations are concerned, there are some 
multidimensional analogs of equation (|8.6(l for the solutions of which the multiple 
elliptic beta integrals on root systems described in this paper may serve as biorthog- 
onality measures. However, their consideration lies beyond the scope of the present 
work. 



Appendix A. Proof of a biorthogonality relation 

We define a pair of functions of z G C as products of two terminating 12V11 
very-well-poised balanced theta hypergeometric series with the argument x = 1 
and different modular parameters: 

Rnm(z) = 12VH — ; — , — , — ,t 3 z, —,q , — ;q,p 

\L4 ^0^4 £]>4 ^2^4 % ^4 

T , fte P P P . h - m Ap" 1 - 1 
>< 12V11 [—;—,—,— ,t 3 z,—,p , — \p,q), (A.l) 

\ 14 to^4 tlM ^2^4 Z C4 

T ^ T/ M*s AAA t 3 _ n Ag™- 1 

\ q t ti t 2 Z U 

( Ab 3 AAA t 3 m Ap" 1 - 1 \ , . . 

\ P to h h z U J 

where n, m e N. Obviously, these functions are symmetric with respect to the 
permutation of p and q. The balancing conditions are used already in these series 
in order to express one of the parameters in terms of the others. 

For m ^ and the fixed parameters q,b r ,r — 0,...,4, the limit as p — > 
is not well defined for R nm {z) and T nm (z). The reason for this comes from the 
quasiperiodicity of 8(z;p) because the limits z — > or z — > oo are not defined for 
it. 

Now, we consider the following integral: 

f dz 
Jmn,kl = / T n i(z)R m k(z)A E (z,t) — , (A. 3) 

Jc mn , kl z 

where A^;(z,t) denotes the weight function l|1.13|l and C mn ,ki is some contour of 
integration. We want to show that a particular choice of C mn ,ki (to be specified 
below) leads to the formula 

Jmn,kl = h n lS mn Skl, (A. 4) 

where h n i are some normalization constants. 

An elliptic extension of the Bailey transformation for the lo^g series, which was 
derived by Frenkel and Turaev in |FT| , has the following form (an alternative proof 
of it was given in |S6|): 

*r u i 4 ^ 0(qto,qs /t4,qso/t 5 ,qto/t i t 5 ;p;q) N 

i2Vu[to;ti,...,tr,q,p) - —77 — —77 7— 7— 

0{qsQ,qt o /t 4 , qt /t 5 , qso/Ut 5 ;p; q) N 

Xi2Vn(so- ) s 1 ,...,s 7 ;q,p), (A.5) 

where l\L=i *m = t 3 Q q 2 , t 6 = q~ N , JVeN, 

qtl soh s t 2 s t 3 , . 

«o = 7-7—, 8i = —, s 2 = — , s 3 = — , (A.6) 
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and S4, S5, se, S7 form an arbitrary permutation of tA,t5,te,t7- Using this identity, 
we can rewrite the functions R m k(z) as follows 

Rmk(z) = n m (p;q)n k (q;p) 

't q q q ± to _ m Aq m ~ 1 



y , 4 q ^ to 

t 4 ti^4 i2^4 ^3^4 z 



;q,p 



where 



x 12V1: 



to. _P P P_ f to 

t 4 t\t 4 t 2 t 4 ^-3^4 z 



Ap 



fc-1 



-;p,q 



(A.7) 



9{qt 3 /t 4 , toh ,t t 2 , A/qt ; p; q) r 



0(qt /ti, ht 3 , t 2 t 3 ,A/qt 3 ; p; q) m ' 
Substituting the explicit series expressions (jA. 2f> and (|A.7|) in the definition of 
Jmn,ki <|A.3|> yields 

n l rn k 



.kl 



y+s' 



where 



K m (p; q)n k (q;p)^2 X! X! X! <7 r+ 'V 

r=0 r'=0 s=0 s'=0 

y d(At 3 q 2r -\tog 2s /U;p) 9(At 3 p 2r ' - 1 ,t p 2s ' /t 4 ; q) 
6(At 3 /q,t /U;p) 9(At 3 /p,t /t 4 ;q) 
B{At 3 /q, A/tp, A/h,A/t 2 , g- n , Aq n ~ x /t 4 ; p; g) r 

6(q, t t 3 , t\t 3 , t 2 t 3 , At 3 q n , hUq 1 -"-^; q) r 
9{At 3 /p, A/tp, A/t^A/h^-^Ap 1 - 1 /^ g;p) r > 

9(p, t t 3 , tit 3 , t 2 t 3 ,At 3 p l , htip 1 - 1 ; q; p) r > 
eltp/U, q/hU, q/t 2 t 4 , q/t 3 t 4 , q' m , Aq m ~ l /t 4 ; p; q) s 
6{q, tohMt^tot^toq^/U, t q 2 - m /A; p; q) s 
^ 9{to/t 4 ,p/t 1 t 4 ,p/t 2 t 4 ,p/t 3 t 4 ,p~ k ,Ap k - 1 /t 4 ;q;p) s , 
6(p, toh,t t 2 , t t 3 , t p k+1 /t 4 , t p 2 - k /A; q; p) s > 

6{zt 3 , z~H 3 ;p; q) r 6(zt , z~H Q ;p; q) ( 



A B (z,t)- 



9(zA, z 1 A;p; q) r 9(zqt 4 1 ,z- 1 qt 4 1 ;p; q) s 
^ 9(zt 3 ,z^ 1 t 3 ;q;p) r : 9(zt , z~Hq- q;p) s > dz 
9(zA,z- 1 A;q;p) r , 9(zpt 4 1 , z- 1 pt 4 1 ;q;p) s , Z ' 
We introduce the notation 



(A.8) 



so that 



to = t q s p s , h=ti, t 2 = t 2 , t 3 = t 3 q r p r 

4 

A=H t m =Aq r p r '. 

m— 

Then, using the transformation property 

9{z;p;q) l = (-zYq l ^/ 2 9(z- 1 q- l + 1 ; p; q) t = 



t 4 = t 4 q s p s 



{-z) l q 



l n l{l-l)/2 



we can rewrite the integral 1A.8|) in the form 



_ ft Q t 4 \ 

J-rs.r's' — 1 I 

pq J 



2ss' 



2rr' ,2(s+s') 
t 4 



q s(s+l) p s'(s' + l) 



C„ 



0{qz 1 ;p;q)-i' 



* / ~x dz 
A E {z,t) — . 

z 



(A.9) 
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But the integral on the right-hand side of (|A.9J) coincides with the elliptic beta 
integral Ijl.l2|l . provided we identify C m nM — T and impose the constraints \t m \ < 
1, \pq\ < \A\. However, the values of the integers r,s,r',s' £ N are not limited; 
starting with their sufficiently large values, we shall have either \A\ — \q r p r A\ < \pq\ 
or |?4| = \q^ s p^ s 't i \ > 1. Now is the moment to specify the contour C m nM- We 
choose it in such a way that formula (|1.12() remains applicable. More precisely, let 
CmnM be a deformation of T such that it separates the poles at z — t i^ 23 p a q b 7 
t 4 p a - k q b ~ m , and A~ 1 p a+1 ~ l q b+1 ~ n , a, b E N, that lie inside C mn M and converge to 
zero, from the poles diverging to infinity, which are obtained from the poles inside 
CmnM by the inversion transformation z — > z . The subscripts m,n,k,l in the 
notation for such a contour indicate that, evidently, the shape of C nm M depends 
on the indices of the functions T„i(z), R m k(z). 
For such a contour C mn ki , we have 



Ir 



, , * 2ss' / , \ 2rr ,2(s+s 
tot 4 \ I ts \ t 4 



pq J \A 



\r 0(^1*3, *2*3, Ht<±\P\ q)r 



9(A/t ,A/t 1 ,A/t 2 ;p; q) r 

x 0(t Q t 3 ;p;q) r+s 9(tot 1 ,t t 2 ,q 1 - r to/A;p;q) s 
9(A/t 4 ;p; q) r+s 9(q/t 1 t i . q/t 2 U, q 1 ^ '/t 3 t 4 ; p; q) s 
9{t 1 t 3 ,t 2 t 3l t 3 t i \q\p) r , 

* ' 6{A/tQ,A/t u A/t 2 )q\p)r> 

9(t t 3 ;q;p) r/+s/ 9(t t 1 ,t t 2 ,p 1 - r ' t /A;q;p) s , 
9(A/t4;q- 1 p) r >+s> 0(p/tiU,p/t 2 t A ,p 1 ' r ' /t 3 ti; q;p) s ' ' 

where the function A/s(t) is fixed in (|1.14() . 

As a result of these manipulations, the quantity J m nM splits into a product 
of two double series each depending only on the indices m, n and k, I separately. 
After an application of the relation (a:p; q) r + s — { a q r \P]q)s{o-\P'iq)r and various 
simplifications, we can write 



where 

J mn (p;q) = K m (p;q)y^q 



JmnM = N E {t)Jmn{p;q)Jki{q\p), 

6{AUq 2r - l ;p) 6(At 3 /q, q~ n , Aq^/U, t 3 U;p; q) r 



9(At 3 /q;p) 9(q,At 3 q n 1 t 3 Uq 1 - n 1 A/U;p;q) r 
ml', (|;^, t t 3 q r , q -^,^-^ 7 q- m ;q,p) . (A.IO) 



The constraint t 2 t 3 — qto imposed upon relation ljA.5l) converts the i2Vn-series 
on its left-hand side into a terminating 10V9 series, whereas on the right-hand side 
only the first term of the corresponding i2Vii-series survives. As a result, we get the 
Frenkel-Turaev sum, or an elliptic generalization of the Jackson sum for terminating 
very-well-poised balanced s^-series: 



10^9(^0! ti, ti, t 5 ,t§, tj; q,p) 

_ 0(qto, gta/tit^, qta/txts, qty/t 4 t 5 ; p; q) N 
9(qt /t 1 t 4 t 5 ,qt /t 1 ,qt /t 4 , qt /t 5 ;p; q) N ' 



(A.ll) 
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where 4^4*5*6*7 = qt% and U = q N , N E N. Application of this sum to the 10V9 
series in (|A. 10|) yields 

y , s Qjqto/h, t 1 t 2 ,At 3 q r - 1 ,q- r ;p; q) rn 

10 9l ' ' '>~ 9{t t 3 , A/qto, Aq r /t 4 , q^/hUip; q) m ' 

Clearly, this expression vanishes for m > r. This means that J mn = for m > n. 

For m < n, we get 

t / \ / , g(At 3 ;p;g) 2 m 9(Aq n - 1 /t 4 ,qt /t 4 ,t 1 t 2 ,q- n ;p;q) m 

J mn [ P ,q) - ^AP,q) e{A/u . p . q)2m 0( hUq l-n jtot ^ A/qtojAt3q n. p . q)m 

x (< 3 t 4 ) m sVriAtaq 2 ™- 1 ; t 3 t 4 , Aq n + m ~ x /t 4 , g™-"; <j,p). 
Applying the summation formula (|A.11J| to the latter gVV series, we get 
K , s 0{Ahq 2m , q m ~ n +\Aq m+n /t 4 , qt 3 t 4 ;p; q) n _ m 

which is equal to zero for n > m due to the factor 9(q m ~ n+1 ;p: q) n -m- As a result, 
Jmnip^q) — h n (p;q)S mn , where the normalization constants have the form 

, , x _ S(A/qt 4 ;p)9(q, qt 3 /U, tph , t t 2l tit 2 , At 3 ;p; q) n q~ n > A ^ 

" W9J ~ 0(A(z 2 ™-Vt4;p)^(lA3i4, to*3,ti*3, *2*3, A/qt 3 , A/qt 4 ; p; q) n 
The fact that J„ m = for ti ^ m provides the desired biorthogonality relation 
QA.4fl . We summarize the result obtained in the form of the theorem that was 
announced in |S4| (it is necessary to apply the series notation introduced in |S5I IS6| 
and permute the parameters t 3 and t 4 in |S4| in order to match with the current 
presentation). 

Theorem 9. Lett m , A E (z,t), J\f E (t) be the same as in Theorem 1. Let C mrit ki de- 
note a positively oriented contour separating the points z = {io,i,2,3P a 9 b j t 4 p a ~ k q b ~ m , 
A~ 1 p a+1 ~ l q b+1 ~ n } a ,beN from the points with the inverse (z — > z^ 1 ) coordinates. 
Then R m k{z) and T n i(z) satisfy the following biorthogonality relation 

f dz 

/ T n i(z)R mk (z)A E (z,t) — = h n iAf E (t)6 mn S k i, (A. 13) 

Jc mn , kl z 

where h n i are the normalization constants, 

9(A/qt 4 ; p)9{q, qt 3 /t 4 , t h, t t 2 , tit 2 , At 3 ; p; q) n q~ n 



h n l — 



9{Aq 2n /qt 4 -p)9{l/t 3 t 4l t t 3 , ht 3 , t 2 t 3 , A/qt 3 , A/qt 4 ;p; q) r 
9(A/pt 4 ; q)0(p,pt 3 /U, t ti,t t 2 , ht 2 , At 3 ; q;p)ip~ l 



" 9(Ap 2l /pt 4 ; q)9(l/t 3 t 4 , t t 3 ,tit 3 ,t 2 t 3 , A/pt 3 , A/pt 4 ; q; p)i' ^ A ' U ' 

As is clear from (|8.5I) and (j8.22jl . we have i? m (z; q,p) — R m0 (z) and T n (z; q,p) = 
T n o(z). These functions R m ,T n are equal to i 2 Vn elliptic hypergeometric series 
with particular choices of the parameters. 

Corollary 10. The functions R m (z\ q,p) and T n (z; q,p) satisfy the following bior- 
thogonality condition 

f dz 

/ T n (z;q,p)R m (z;q,p)A E (z,t) — = hJ\f E (t)6 mn , (A.15) 

Jc„,„ z 



where the constants h n are fixed in fA.lty) and the contour C mn encircles the poles 
of the integrand located at z — {io,i,2,3<Z a ?A t 4 p a q b ~ m , A~ 1 p a+1 q b+1 ~ n } a ^eN and 
separates them from the poles with inverse z — > z^ 1 coordinates. 
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The biorthogonal rational functions R m (z;q,p) and T n (z;q,p) describe elliptic 
generalizations of the Rahman set of continuous lo'&g functions |R1| to which they 
are reduced in the limit as p — > 0. Accordingly, in this limit, formula (|A.15|I is 
reduced to the Rahman biorthogonality condition. 



Appendix B. Integral representations for 12-E11 series 

Here we derive an integral representation for the product of two terminating 
12-E11 (more precisely, 12 Vu) series with some particular choice of parameters. For 
this, we apply an elliptic generalization of the technique used in |R,1| for the deriva- 
tion of the contour integral representation for a terminating io<l>9 series. 

Theorem 11. Suppose that five parameters tk,k — 0, ...,4, satisfy the condi- 
tions of Theorem 1. Denote by m,n two positive integers and by C mn a positively 
oriented contour such that for all a,b £ N it separates the points z — {tf ! p a q b , 
j^-lqb+l— mpa+l-ra| j rom iheir partners with the inverse coordinates (z — > z^ 1 ). 
Under these conditions, the following integral representation for the product of two 
12V11 terminating very-well-poised balanced theta hypergeometric series at X = 1 
holds true: 

(At _ m A 2 q m - 1 

\iv\\ — ; ot, ton, <o^2, ^o*3; toHi q > \q,p 

\ q a 

X12V11 ( — ;{3, t a t 1 ,t t 2 ,t t 3l tQt 4l p' n , ;p,q 

_ 1 9(At ,^;p;q) m 9(At ,±;q;p) n 

A ( M^^ z ;p;<iUe(f,£;q;p)n dz 

L " 9(Az±;p;q) m 9(Az±;q;p) n z ' 



where a and [3 are arbitrary complex parameters. 

Proof. Under the conditions imposed upon the parameters in the formulation of 
this theorem, the following relations are true: 



/ A £ (z,t) 



A 



9(zt ,z 1 t ;p;q) l 9(zt ,z 1 t ;q;p) j dz 
9(zA, z- l A;p\ q) l 9(zA, z^A; q;p)j z 

Af E (toq*p> ,ti, ...,t A ) 



9(t t!, . . . ,toU;p;q)i9(tot 1: . . . ,t t 4 ;q;p)j r , . , , 

9(A/t u A/U-p- q) l 9(A/t 1 , A/U; q;p) j E{ >' l ' J 



We multiply (|B.2|) by the factor 

^{AW^-p) 9(At /q, a, <T™ A^^/a-p; q) t 



xp 



9{At /q-p) 9{q, Ato/a, At Q q m , aqi-^to/A; p; q) t 
j 9(At oP ^- 1 ; q) 9(At /p, (3, p~ n , A 2 p n ~ 1 /P; g; p) 3 



9(At /p; q) 9{ P , Ato/0, At oP n , p P *-n /A; q-p) 3 
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where a and (5 are arbitrary complex parameters, and sum over i from to m and 
over j from to n. As a result, we get the relation 

A E (z, t) 10 V g (At /q; toZ, hz- 1 , a, q~ m , A 2 q m ~ 1 /a; q,p) 

dz 

x ioV9(Ato/p;t z,toz~ 1 ,l3,p~ n ,A 2 p n - 1 /a;p,q) — 

= i2Vn(At Q /q;a,q~ m ,A 2 q m - 1 /a,tot 1 ,...,t t 4 ;q,p) 

x i2V n (At /p; f3,p~ n , A 2 p n - 1 /I3, t h, t Q U; p, q)M E (t). 

Application of the Frenkel-Turaev sum to the 10V9 series standing under the integral 
sign leads to (|B.1|) . □ 

For n — we get an integral representation of a single terminating 12 V\\ series, 
which can be reduced further to the 10 $9 g-series level |Rlj by letting p — » 0. 
However, for n ^ the limit as p — > is not well defined, and formula IjB.lf) exists 
only at the elliptic level. 
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